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Abstract — Representing the relativistic physical fields as sections of the Clifford Bundle (or of the 
Spin-Clifford Bundle) of Minkowski spacetime we show that all the relativistic wave equations satisfied 
by these fields possess solutions traveling with arbitrary speeds < V < OO. By giving rigorous 
mathematical definitions of reference frames and of the Principle of Relativity (PR) we prove that 
physical realizations of the V > 1 solutions of, e.g., the Maxwell equations imply in a breakdown of 
the PR, but in no contradiction at all with known physical facts. 



1. INTRODUCTION 

In this paper we present methods for constructing solutions with arbitrary speeds 
(0 < v < oo)Q of the main relativistic wave equations of physics, namely the (scalar) 
homogeneous wave equation (HWE) and Maxwell, Weyl, Klein-Gordon and Dirac equa- 
tions. Several examples are worked in detail. 

In Section ^| we show how to represent all fields mentioned above as sections of 
the Cliford bundle C£(M) or of the Spin-Clifford bundle Ct Sv in + (1,3) (M) of Minkowski 
spacetime jl], |[ |). We then show that if the HWE and the Klein-Gordon equation 
(KGE) have solutions with arbitrary speeds < v < oo, it then follows that also 
Maxwell, Weyl and Dirac equations have such solutions. 

The solutions we are going to exhibit for the relativistic wave equation are families of 
undistorted progressive waves (UPWs) . By UP W, following Courant and Hilbert Q we 
mean that the waves are distortion-free, i.e. that they are translationally invariant and 
thus do not spread, or that they reconstruct their original form after a certain period 
of time. 

In Section || we present < v < oo UPW solutions of the HWE and of the KGE. The 
meaning of the velocity of propagation of a UPW and the concepts of phase and group 
velocities are discussed in details. Misconceptions regarding the velocity of propagation 
of energy of a given wave are clarified. In Section |J we present w ^ 1 UPW solutions 
of Maxwell equations, discussing the remarkable characteristics of these solutions in 
contrast to the v = 1 solutions. One of the principal novelties is that the v ^ 1 

1 We use units such that c = 1, where c is the so-called speed of light in vacuum. 
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solutions have in general non-null field invariants and are not transverse waves. We also 
discuss in this section the important problem of the velocity of transport of energy for 
the v > 1 UPWs solutions of Maxwell equations. 

It must be said that there are experimental data |^|, |6| obtained with techniques 
developed by Lu and Greenleaf |?], ||] showing the existence of pressure waves moving 
with speeds^] v < c s and v > c s which confirm the theoretical predictions, showing in 
particular that the energy associated with the pressure waves can travel with speed 
v £ > c s . In || it is discussed the possibility of designing physical devices for launching 
in physical space the v ^ 1 solutions of Maxwell equations. 

If this is indeed possible we must investigate the status of the Principle of Relativity 
(PR). This is done in Section ||, where we give a thoughtful and rigorous mathematical 
definition of this principle and its relation with Lorentz invariance. We prove that the 
v > 1 solutions of Maxwell equations imply necessarily in a breakdown of the PR, 
which is necessary in order to avoid logical contradictions. In Section || we present our 
conclusions. 



2. A UNIFIED THEORY FOR CONSTRUCTION OF UPW SOLUTIONS OF MAXWELL, 
DIRAC AND WEYL EQUATIONS 

To fix the notations we recall here the main results concerning the theory of Clifford 
algebras (and bundles) and their relationship with the Grassmann algebras (and bun- 
dles) . In particular a self-consistent presentation of the so called spacetime and Pauli 
algebras is given. Also the concept of Dirac-Hestenes spinors and their relationship with 
the usual Dirac spinors used by physicists is clarified. We introduce moreover the con- 
cept of the Clifford bundle of spacetime and the Clifford calculus. As we shall see, this 
formalism provides a unified theory for the construction of UPW subluminal, luminal 
and superluminal solutions of Maxwell, Dirac and Weyl equations once we have arbi- 
trary speed (0 < v < oo) solutions of the homogeneous wave equation (HWE) and of 
the Klein-Gordon equation KGE. More details on these topics can be found in @, §,§). 

2.1. Exterior, Grassmann and Clifford algebras 

Let V be a n-dimensional real^j vector space, V* its dual space, T r V the space of r- 
contravariant tensors over V(r > 0,T°V = M, T X V = V) and let TV be the tensor 
algebra of V. 

We recall that the exterior algebra of V is the quotient algebra 



/\V = TV/J (2.1) 

where J is the bilateral ideal in TV generated by elements of the form u ® v + 
v <£> u,u, v £ V. The elements of f\V will be called multivectors, or multiforms if 
they are elements oi /\V* . 

Let p : TV — > /\V be the canonical projection of TV onto /\V. Multiplication in /\V 
will be denoted as usual by A : f\V — ► /\V and called exterior (or wedge or Grassmann) 
product. We have 

AAB = p(A ® B). (2.2) 



2 Here c s is the speed of sound in water. 

3 Here IR denotes the real field. 
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We recall that f\V is a 2"-dimensional associative algebra with unity.^j In addition it 
is a Z-graded algebra, i.e., 

/\V = ®A r V; A r V A A S V c A r+s V ; (2.3) 

r, s > 0, where [\^V = p(T r V) is the (")- dimensional subspace of r- vectors, /\ V = M, 

A V = V; = {0} if r > m - K A G f° r some r(r = 0, . . . , n) then A is said 
to be homogeneous, otherwise it is said to be inhomogeneous. For A p £ A P V and 
B q e A 9v we have 

A p AB q = (—l) m B q A A p . (2.4) 
Let (Ei, . . . , -E„) be a basis for y. Then a basis for A^ is 
{1, Ei, ...E n , Ei AE 2 , Ei AE n , ... , En-iSE ni Ex A E 2 AE 3 , . . . , Ei AE 2 A . . . A E n } (2.5) 
Then, if A € A^j we can write 

A = S + A l E, + j { A ij Ei A Ej + ^ A ljk E l AEjAE k + ... + PE X A E 2 A . . . A E n , (2.6) 

where 5, A ij ,A^ k , . . . , P e M and A ij = -A ji , etc... The element 

E n+ i =E 1 AE 2 ...AE n (2.7) 

is called the pseudoscalar of the algebra f\V . (The analogous element for A^* i s a i so 
called the volume element). We define the projector Qf. : /\V — > A ^ by 3^ H 

n 

A fc ,for A = £\4 Jb ,4 b G /\ ky - 

Now let g E T 2 V* be a metric for V of signature (p,q), i.e., g : V x V — > 1R and 
let g^ 1 G T 2 V,g~ 1 e : V* x V* — > _ZR be the metric of the dual space. If it, v G V and 
a,(3 eV* such that a (it) = l,/3(u) = 1 we have 

g(u,v)=g-\a,P) . (2.8) 

We can use g to induce a scalar product on A^> G : A^ x — * We dchnc 

G(A,B) — det(g(ui,Vj)) (2.9) 

for homogeneous multivectors A = ui A . . . A u r € A^i B = ui A . . . A u r € A^j 
iti,Uj € A^j *i J ' = 1) ■ • ■ > r - This scalar product is extended to all V due to linearity 
and associativity. G(A, B) = if A e A'X B € A*^> r 7^ s - When both A, B e f\V, 
G(A, B) means the product AB. The algebra A^ endowed with this scalar product is 
called Grassmann algebra and will be denoted A {Vid)- 

On A^ and A(^iff) there are two important involutivc morphisms: 

(i) Main automorphism A : f\V — ► A^i 

(A A By = AAB, A,Be AV ; 
A — A if AeA°V, A = -A if AetfV. ' 



4 A V is what old physics textbooks call the algebra of antisymmetric tensors. 
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(ii) Reversion ~ : /\V -> /\V, 

(A A B)~ = BAA, A,Be f\V ; 



A = A if ieA'VffiA 1 ^- (2 ' n) 

We define also: 

(iii) Conjugation: — : /\F — > /\V, 

3= = (1) A , VA e /\V. (2.12) 

We introduce now the important concepts of left contraction J : f\V x f\V — > /\V 
and right contraction |_: A^ x ~~ * A^ through the definitions 

G(A\B,C) =G(B,AAC); G(A[B,C) = G(A,C A B), VC e A^ ■ (2-13) 
J and [ satisfy the rules 

1. x\y = x.y = g(x,y)\ x[y = x.y = g(x,y) ; 

2. x\(AAB) = {x\A)AB + AA{x\B)- 

3. {AAB)\x = A A (B[x) + (A[x) AB ; l ' ' 

4. (AAB)JC = AJ(BJC); AL(-B AC) = {A[B)[C ; 

where /\V, i,B,Ce A^- 

The notation A_B will be used for contractions when it is clear from the context 
which factor is the contractor and which factor is being contracted. When just one of 
the factors is homogeneous, it is understood to be the contractor. When both factors 
are homogeneous we agree that the one with the lowest degree is the contractor, so that 
for A £ /\ r V, B e /\ S V we have A.B = A\B if r < s, A.B = A[B, if r > s. From the 



definitions and eq.(2.14) we easily verify that 

s\A = 0, A[s = Vs e M, A e /\V . (2.15) 

We are now ready to present the definition of the real Clifford algebra C£(V, g) asso- 
ciated with the pair (V,g). In order to do that we define the Clifford product (denoted 
by juxtaposition of symbols) between x E V and A € A^ by 

xA = x\A + xAA = x.A + xAA 

and extend this product by linearity and associativity to all of A^- 

Equipped with the Clifford product A^ becomes isomorphic to the Clifford algebra 
C^(V,.g)J^] Observe that A^ equipped with the exterior product and C£(V,g) equipped 
with the Clifford product are, of course, not isomorphic as algebras. However, A^ and 
C£(V, g) are isomorphic as linear spaces over 1R. 



Consider the basis of A^ given by (2J5) and suppose that 



~T~1 , 1 — J = 1, , J), 

<!l£ i .Ej) = { -1, i = j=p+l,. , ..,p + q; (2.16) 
0, otherwise. 



5 We can show that C£(V,g) = TV/ J where J is the bilateral ideal on TV generated by elements of 
the form a®b+b®a- 2g(a, b), a, b e V C TV jlj. 
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Then it is clear that EiEj = Ei.Ej + Ei A E 3 — Ei A Ej for all i ^ j. Since f\V and 
C^(V, 5) are isomorphic as linear spaces we can write for X £ C£(V, g) 

X = S + XiE* + -X ij E i E i + j^X ljk E i E 3 E k + ■■■ PE 1 E 2 E 3 ...E n , (2.17) 

where E % .Ej — 6* and {E 1 }, i — l,...,n is called the reciprocal basis of V. Also 
S, Xi, Xij , . . . , P £ Ft and Xij = —Xji, etc. 

For g of signature (p, q) as in eq.( 2.16| ) C£(V, g) is denoted C£ p>q . Using the projector 



operator ()& defined above we can show that the contraction A r .B s , for A r £ f\ V C 
C£ p<q , B s £ f\ s V C C£ Piq is given by 

A r .B s = { i A f°^-s\ if r,s>0 

' if r = or s = v ; 



Eq.(2.18) defines then an inner product in C£ 



We now define the Hodge star operator * : CtiV, g) — » C£ (V, g) by 

*A = AE n+l , A £ C£(V,g). (2.19) 

A simple calculation shows that *|yy p maps f\ p V — > f\ n p V for p = 0, 1, . . . , n. We 
observe that C£ PA is a ^-graded algebra. This means the following. Let C£ pq {C£~ q ) de- 
note the set of even (odd) multivectors of C£ Piq , i.e., elements of f\ 2r V c C£ p , q (/\ 2r+1 c 
C£ p , q )- We have C£ pq C£ pq c C£ pq , C£ pq Cl pq c C£ pq , C£ p q C£ p q c C£ pq , C£ p q C£ p q c 
C^~ g . The C£+ ? is a sub-algebra of C£ Pi(? , called the even sub-algebra of Ci? Pjg . All 
Clifford algebras C£ Ptq are isomorphic to matrix algebras over the fields M, <S or M , 
respectively the real, complex and quaternion fields. We can find in [^j tables giving 
the representations of C£ Ptq as matrix algebras. For what follows we need to know the 
following representations: 

Complex numbers — C£q.i ~ <S 
Quaternions — C£q,2 — M 

Pauli algebra - C£ 3fi ^ M 2 (d7) 



Spacetime algebra — C£\$ ~ M 2 (1H) 
Majorana algebra — C£s t \ ~ M^M) 
Dirac algebra - C4,i — M 4 (d7) 



(2.20) 



Since it is a theorem that C£~£ q ~ C£ q ^\ for p > 1 and C£+ 9 ~ C£ p q ^\ for g > 1 we 
have the following useful identifications to be used later: 

«+ 3 ~ c# tl ~ C4,o; ~ C4, 3 ; «+ ~ h ; H+ ~ <s. (2.21) 

A very important result is that the Dirac algebra is the tensor complexification of both 
C£i t 3 and C£a t i, i.e., 

Ch,i ~ d7 ® «i, 3 ; « 4 ,i ~ (C* ® C4,i (2.22) 

Since it is a well known result that 1H is represented by a subset of invertible two by 
two complex matrices belonging to M2((F), eq.(2.21) and eq.( [2.22j ) show that C£i t 3 has 



also a complex 4x4 matrix representation which can be made identical to the M${(E) 
representation of C£i \ . 
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Before ending this section we recall that the Clifford product between two general 
elements A, B S C£ Ptq can be written 

AB = ^{A) r {B) s = ^A r B s 

r,s r,s 

= ^2((ArB a ) lr _ sl + (A r B s ) lr _ sl+2 + ■■■ + (A r B s ) r+s ). (2.23) 

We define also the norm of a multivector A E C£ p ^ q by 

\A\ 2 = (AA) . (2.24) 
If A £ Ci^fl is homogeneous and \A\ 2 ^ 0, the inverse of A is^] 

A- l = A/\A\ 2 , A- 1 A = AA- 1 = 1. (2.25) 
2.2. The spacetime and Pauli algebras 

We call SI 1 ' 3 = (R 4 ,g) where g : R 4 x R 4 -> R is a Lorentzian metric of signature 
(1, 3). R 13 is called the Minkowski vector space. Let {E^}, fj, = 0, 1, 2, 3 be a basis of 
R 4 ; we have 

f +1 (i = v = 0; 
= 7?^ = ^ -1 /x = i/=l,2,3; (2.26) 
[ otherwise. 

The fundamental relation generating the spacetime algebra Cl\^ is then 

EpE v +E v Ep = 2Ti lu ,. (2.27) 



Eq.(2.27) is identical to the relation satisfied by the famous Dirac (gamma) matrices 
and indeed we know from Section 2.1 that the E^ have a complex matrix representation 
in Mi( <&). Naturally, dimCfi^ = 16. The pseudoscalar of Cl\^ will be denoted 
by E$ = E0E1E2E3 and E 2 = — 1. E§ anticommutes with odd multivectors and 
commutes with even multivectors. We call {E 11 }, ri = 0,1,2,3, such that E^ l .E v = 5£ 
the reciprocal basis to {E^}. 

We call the pair R 1 ' 3 * = {R 4 ,g- 4 ) the dual space of (R 4 ,g) and call {1^}, 11 = 
0, 1, 2, 3 the dual basis to {E^}. Analogously, {T^}, fj, = 0, 1, 2,3 such that g^ x {T^, T v ) = 
8^ is the reciprocal basis to -jT^}. The Clifford algebra associated to R 1 ' 3 * will be de- 
noted *C£\ t 3 ~ Cii^ ~ M2(JH). Of course we have the fundamental relation 

T^T" + rT^ = 2rf", (2.28) 

where rf' v = r/^. 

The Pauli algebra C^o is the Clifford algebra of (R 3 , <?b), i.e. of Euclidean space 
equipped with the Euclidean metric gE- If {^i}, i = 1,2,3 is an orthonormal basis of 
R 3 , i.e. £j.£j = gE(^i,^j) = then the Clifford algebra Cl^^ is generated by the 
fundamental relation 

EiEj +£,,•£; = 2%. (2.29) 



3 The calculation of A 1 (when it exists) for a general A £ C£ p , q is not so simple. 



On UPW solutions of relativistic wave equations 



361 



/ = EiE 2 E 3 is the pseudoscalar of CI^q. We verify that J 2 = —1 and that / commutes 
with all X e C^3,o, so that I is like i = \/^T. A basis for C£ 3i q is (1, Ej, EiE,-, EiE 2 E 3 ). 
Taking into account that 

* (EiEa) - E 2 Ei7 - E 3 ; *(Ei£ 3 ) = E3E1/ = -E 2 ; *(£ 2 £ 3 ) = E 3 E 2 7 = E i; (2.30) 

we see that VXe C^3,o can be written as 

X = (s + Ip) + (A + IB) s,p € M; A,B e R 3 ; (2.31) 

i.e., X i s "formally the sum" of a "complex number" and a "complex vector". From 
eq.( 2.21 ) we see that Ciz.o — C£^ 3 . We can exhibit this isomorphism by identifying 
Efe = Ek.Ev — Efc, k — 1,2,3 and where E 2 = 1 is timelike. We define in C£ 3 ,o the 
operator of spatial inversion. For X € Cl^fi as in eq.( [2.31 ), 



X i-> X* = (s + Ip) -{A + IB). 



(2.32) 



With the above identification p — p^E^ £ /\ (R 1 ' 3 ) C Cl\ 3 can be represented in 
«+ 3 ~« 3 ,oby 

p^pE Q =p° + p% =p° +p. 



For / = y^Ef.Eu e /\ a (M l > 3 ) C Cii,s where 



we can write 



f 


-£1 


-£2 


-£s \ 







-S3 


B 2 


£2 


S3 





-Bi 




-62 




/ 


EiEj 


= -/° 


% + 


/* J E,Ej 



(2.33) 



(2.34) 



and taking into account ( |2.29 ) and ( 2.30| ) we can write 

f = £ + IB / = EiE 2 E 3 . 



(2.35) 



Eq.( 2.35 ) shows that any bivector F £ /\ (R 1 ' 3 ) C C^i, 3 can be represented in Ci? 3j o — 
C£f 3 by a sum of a vector and a Pauli bivector or a "complex vector" . 

We define next the Spin groups Spin + (3,0) ~ SU(2) and Spin + (1,3) ~ SL(2, <D), 
which are respectively the covering groups of SO '+(3), the special rotation group and 
SO+(l,3) ~ £l_, the restricted orthochronous Lorentz group. We have 



Spin + (3,0) = { J Re«+ | |i?| = l}, 
Spin + (1,3)={C/G^+| \U\ = 1}. 



(2.36) 
(2.37) 



An arbitrary Lorentz rotation is given for a G Cl\$- by a ^ UaU — UaU , 
U € Spin, (1,3) We can prove that any U € Spin + (1,3) can be written in the form 
U = ±c^, / e /\ 2 (R}' 3 ), and the choice of the sign can always be positive except when 
U = —e* with f 2 = 0. When f 2 > 0, [7 is a boost and when / 2 < 0, {/ is a spatial 
rotation. We end this Section with the definitions of minimal left ideals of C£ Ptq and of 
geometrically equivalent ideals. 
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We say that e G C£ PtQ is idempotent if e 2 = e; it is called a primitive idempotent if it 
cannot be written as a sum of two mutually annihilating idempotents, i.e., e^e' + e", 
with (e') 2 = e', (e") 2 = e", e'e" = e'V = 0. 

The sub-algebra I e C C£ Pi9 is called a Ze/t idea? of C£ Ptq if V"0 € I e and VX € C£ Pjg we 
have Xtp G 7 e (a similar definition exists for right ideals). 

An ideal is said to be minimal if it contains only trivial sub-ideals. It can be shown 
that the minimal left ideals of C£ Piq are of the form Ci p , q e, where e is a primitive 
idempotent. 

Consider now C£\ t z and the orthonormal bases S = {E^} and X = {E^} where 
= UE^U, U G Spin + (1,3). We can easily verify that the following elements are 
primitive idempotents of Cfi 3: 



(2.38) 



e E = i(l + £ ); 4 = i(l + £ 3 £o); 4 = £(l + £i3»£fe); 
e ± = ±(l + £ ); e^ = ±(l + £ 3 £ ); = §(1 

It is trivial to verify that e E and e^. are related by 

e ± = UesU' 1 , u G Spin + (1, 3). (2.39) 

There is no element U E Spin + (1, 3) relating, e.g. e' s with ejj. Consider now the ideals 
Is = C£i ; 3es and ij, = Cli^e y;. W e say that Is and 1^ are geometrically equivalent 
if es and ej, are related by eq.( 2.39| ). Since 1% — C^i )3 es and ij, = Ci\^e-^ and since 
C4, 3 [/ ~ «i, 3 W G Spin + (1, 3), we can write 

I ± = IvU- 1 (2.40) 



Eq.( 2.40 ) defines a correspondence between elements of ideals that are geometrically 
equ ivale nt. The quotient set {/s}/7£ where 1Z is the equivalence relation given by 
eq.( 2.40| ) is called the space of the Dirac algebraic spinors S. Of course in the basis £ 
the spinor is represented by ip s and 

= ^U- 1 (2.41) 



Section 2.4, where we introduce the fundamental concept of Dirac-Hestenes spinors, 
will clarify the meaning of the above definitions. 

2.3. Dirac algebra Ci^i, its relation with C£± t s and Dirac-Hestenes spinors 

Consider the vector space M 4 ' 1 = (M 5 ,g) and let {E a }, a = 0, 1, . . . , 4 be an orthonor- 
mal basis: 

f 1, = 6=1,2,3,4; 
g{E a ,E b ) =1 -1, a -6 = 0; (2.42) 
[ 0, otherwise. 

Let C£i t i be the Clifford algebra of Ft 4 ' 1 and I = E0E1E2E3E4 the corresponding 
pseudoscalar. Note that I 2 — — 1 and that E a I = IE a , and thus M(B IM is the center 
of C£^i and the pseudoscalar / plays therefore the role of the imaginary unit (as I in 
the case of the Pauli algebra) . 
Let us define 

r M = £ M £ 4 , M = 0,l,2,3. (2.43) 

Then, 

r^r, + = 2^. (2.44) 
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One can easily see from (2.43), (2.44) and with I playing the role of the imaginary unity 
(i = \J — 1) that Cl^i is isomorphic to the complexified spacetime algebra, i.e., 



C£ 4 ,i ^ <S®C£ h3 . 



(2.45) 



Indeed, each X € Ctn can be written 



X 



a + X a E n 



1 



f 



1 



TP TP i vabc TP TP TP 

= (a + 76) + JCE^, - (Y^ 1 Efj,I) - IX S E E 1 E 2 E 3 

1 



X abcd E a E b E c E d + lb 



(a + 76) + (A" + ib")^ + -(A" v + jb"")^^ + • • ■ + P r r 1 r 2 r 3 



(2.46) 



a,b,X^,Y^,...,B^,... G J?. 

Moreover, the even sub-algebra of Ci?4,i is isomorphic to dig, i.e., 

C^4_i, the complexified spacetime algebra is the well known Dirac algebra studied in 
physics textbooks. Indeed Cl^.\ is isomorphic to M^dJ), the algebra of 4 x 4 matrices 
over the complex. One representation (the standard one) of the Y ^ defined by eq.( 2.43j ) 
is 



/ 1 \ 

10 

0-10 

\ -1 ) 

I i \ 

-i 

-i 

\ i J 



(2.47) 



/ -1 \ 
0-10 
10 

y i o o o / 

/ -1 \ 
1 
10 

y o -l o o / 

TTi^) = e s i(l + 7T 1 r 2 ) where 

es = — (1 + To) is a primitive idempotent of Cl\ x ~ Cl\$. It generates the left minimal 
ideal I± 1 = Ci^xf and we can easily verify by explicit computation that 



Consider the idempotent / 
1 



i(i + r )i(i 



(2- 



Consider now the usual Dirac spinor \ty) 6 <S A . There is an obvious isomorphism 
between (F 4 and minimal left ideals of M^fS), given by 



(F 4 9 |*) 



V ^ j 



( i/j x o o o \ 

i> 2 o o o 

V> 3 o o o 

\ipi J 



* S minimal left ideal of M 4 (d7) 

(2.49) 
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One can, of course, work with <3> instead of and since M^tE) is a representation 
of the Dirac algebra Ct^i ~ (V ® C^i,3 we can work with the corresponding ideal lf 1 of 
the Dirac algebra. The isomorphisms discussed above tell us that 

III = C£ iA f = (d7 <8> C4 >3 )/ - C<i/ ~ C*i, 3 / = («i, 3 e E )i(l + JTxra). (2.50) 

Note that in the last equality we have a minimal left ideal Cii^e-s of the spacetime 
algebra. Moreover we have 

// = r 2 r x /. (2.5i) 

These results mean that we can work with the ideal Clx^e-Y, once we identify r 2 Ti as 
playing in C£\ t3 the role of the imaginary unit. We can verify by explicit calculation 
that 

di, 3 e s =Ce+ 3 e s . (2.52) 

We see that what the idempotent makes is to "kill" redundant degrees of freedom. Since 
dimC-^3 = dimC-^3,0 = 8 we can work with Cl\ 3 instead of C^i )3 es (this is not the case 
for Cl\fi or <£ ® C4,3 since dimC^i, 3 = 16 and dim( <£ (g> C£i }3 ) = 32). We have thus 
established the isomorphism 

d7 4 ~ C£+ 3 . (2.53) 
We shall call ^ s the representative of |*) in Cli \. It is related to ife £ {Cl\ J + ~ Cl\ z 

by 

*E = VE^(i + r )^(i + ir 1 r 2 ). (2.54) 

Such a ips will be called the "representative" of a Dirac-Hestenes spinor in the basis 
(r ,ri,r 2 ,r 3 ) of Cl\ x - C£i, 3 - Its standard matrix representation is 



4>i 



( V>i 




V>3 


V>4* \ 


1p2 




V>4 


r 3 


^3 




Vl 




V V>4 




V>2 


n J 



(2.55) 



A Dirac-Hestenes spinor is an element of the quotient set C£i i3 /1Z such that given 
two orthonormal basis E = {1^}, E = {1^} of K 1 ' 3 C C£i, 3 , ipv £ C^ 3 , then ~ 
^(modft) if and only if iP t = ipsU' 1 with E = £(E) = [/EC/" 1 , [/ € Spin + (1,3), 
£ S SO+(l,3) and W(J7) = C where W : Spin + (1,3) -> 50+ (1,3) is the universal 
double covering of 5*0+ (1,3). We already said that ips is the representative of the 
Dirac-Hestenes spinors in the basis E. When no confusion arises we shall write only 
1 



ip instead of ips. ips — (1 + To) is the Dirac algebraic spinor introduced in Section 2.2 



From now on we work with C£\_ 3 . Then ip £ Ct\ 3 can be written as 

ip = s + f + r 5 p, r 5 = r r 1 r 2 r 3 , (2.56) 

S,P eM, and F e A 2 (R 1 ' 3 ) C C£ h3 is a bivector. 

Suppose now that ip is nonsingular, i.e., ipip ^ 0. Since ip £ C£f 3 we have 

ipip — (T -\- r 5 w, a,u eM. (2.57) 



Define p = \/ a 2 + w 2 , tan/3 = w/cr. Then we have 

ip = y/pf'WR, (2.58) 
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where R £ Spin + (1, 3), p £ JR + and < f3 < 2ir is the so called Ivon-Takabayasi angle. 
This is the canonical decomposition of Dirac-Hestenes spinors and reveals the secret 
geometrical meaning of spinors, for if X £ M 1,3 C C£\$ 

iPX^ = pRXR = pY , Y = RXR £ M 1 ' 3 C C£x, 3 , (2.59) 

i.e., a Dirac Hestenes spinor acting on a vector produces a Lorentz rotation plus a 
dilation of the vector. 

A Weyl spinor ip £ Cli^/H is such that its representative in a given frame £ satisfy 
the condition [[To) 

j 5 ip = ±^721 ■ (2.60) 

Such spinors are called positive and negative eigenstates of 75 and are denoted by ip+ 
(ip-)- For a general ip £ (%f 3 we can write 

^± = 2^^75^721]- (2.61) 

2.4. The Clifford bundle of differential forms and the Spin-Clifford bundle 

Let M. = (M, g, D) be Minkowski spacetime, where (M, g) is a four dimensional 
time oriented and spacetime oriented Lorentzian manifold, with M ~ _ZR 4 and with 
g £ sec(T*M x T*M) being a Lorentzian metric of signature (1,3){] T*M[TM] is 
the cotangent [tangent] bundle. T*M = U xBM T*M [TM = U xf z M T x M] and T X M ~ 
T*M ~ M 1 ' 3 , the Minkowski vector space already defined above. D is the Levi-Civita 
connection of g, i.e., D(g) = 0, T(D) = 0. Also = 0, T and K being re- 

spectively the torsion and curvature tensors. Now, the Clifford bundle of differential 
forms C£(M) is the vector bundle of algebras Cl(M) = U xeM C£(T*M) where Vx £ M, 
C£(T*M) ~ C^i^, the spacetime algebra. As a linear space C£(T*M) is isomorphic to 
the exterior algebra /\(T*M) ~ /{(M 1 ' 3 *) of the space ffi 1 ' 3 * dual of JR 1 ' 3 . Then the so 
called Cartan bundle /\(M) = \J xeM /\{T*M) can be thought as "embedded" in C£{M). 
In this way sections of C£(M) can be represented as a sum of inhomogeneous differential 
forms iJ. 

Let {e^}, p = 0, 1, 2, 3, e M € sec TM be an orthonormal basis of TM, i.e., g{e^, e„) = 
tyiv = diag(l, —1,-1,-1) and let its dual basis be {7^}, \i = 0, 1, 2, 3, 7^ £ sec /\ (M) C 
secC£(M). Then if g" 1 £ sec(TAf x TM) is the metric on T*M, we have g" 1 ^, 7") = 
rf- v = diag(l, —1,-1,-1). The fundamental Clifford product is generated by 

yy + yy* = 2^*". (2.62) 

We introduce also the reciprocal basis {7^}, p = 0,1, 2, 3, 7^.7" = 5", 7^ £ sec /\ (M) C 
sec C£(M). Then C £ sec C£(M) can be written as 

C = s + + i&"*Vr„ + ^a^'wyp + P75, (2.63) 

where 75 = 70717273 and v, , a» vp , p £ sec/\°(M) C secC£(M). 

Besides C£(M) we need to introduce another vector bundle, Cfgpin, (i,3)(M), called 
the Spin-Clifford bundle j^, [|, which is a quotient bundle,?. e., C£s pin ,(1,3) = C£(M)/H. 
This means that the sections of C^spin (1,3) (Af) are some special equivalence classes of 

7 Here sec means Section of a given bundle. 
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sections of the Clifford bundle, i.e., they are equivalence sections of non-homogeneous 
differential forms. A given Section if> G secC£g p i n+ (i 3 ) is then represented by ips, 

. . . G secC£(M) where E, £, . . . are orthonormal bases of /\ (M) C C£(M), with 

yj t = ( 2 - 64 ) 

and Vx G M, R{x) G Spin + (1, 3). 

Dirac-Hestenes spinor fields (DHSF) are sections of C£^ pin ^ 13 ^(M), the even sub- 
bundle of C€ Sp in + (i,3)(^)- The representative of a DHSF on C£(M) in the basis £ is 
then 

V>£ = s + + P75 (2.65) 

The Hodge star map * : f\ p (M) — ► f\~ p (M) can be represented in C£(M) by the 
following algebraic operation: 

*A = A l5 , (2.66) 

for A G aec/\ p (M) C «(M). 

Let d and <5 be respectively the differential and Hodge codifferential operators acting 
on sections of f\(M) C Ci{M). We have 

d:AW-A P+1 (M)- 

If sec A P ( M ) 3 w p = —,u a R . . . 7 Q A 7^ A . . . then, 
p! 



dw p = ie,K„ . . .)t" A 7 Q A 7^ A . . . (2.67) 

and d 2 = 0. Also, 



8 : /\ P (M) -> /\ p (M), du p = (-) p *- 1 d*u p , (2.68) 

where = = identity and 5 2 = 0. 

The Dirac operator acting on sections of C£(M) is the invariant first order differential 
operator d : C£(M) -> «(M) 

9 = ^D e , (2.69) 
and it holds the very important result (see e.g. 

<9 = <9A+<9. = d-<5 (2.70) 

For Lu p G sec f\ p {M) C sec C£(M) 



The operator □ = 99 = (d — <5)(d — <5) = — (d<5 + <5d) is called Hodge Laplacian. 



(2.71) 
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2.5. Maxwell theory in Ct(M) and the Hertz potential 



We shall need the concepts of incrtial reference frames (I), observers and naturally 
adapted coordinate systems. 

Let M = (M, g, D) be Minkowski spacetime. An inertial reference frame (irf) / is 
a timelike vector field I £ sec TM pointing into the future such that g(I,T) = 1 and 
DI = 0. Each integral line of I is called an inertial observer. The coordinate functions 
(x^), fj, = 0, 1, 2, 3 of a chart of the maximal atlas of M are said to be a naturally 
adapted coordinate system to I (nacs/J) if I = d/dx° pi Q. Putting / = eo we can 
find ei = d/dx\i — 1,2,3 such that g(e M ,e„) = i]^ and the coordinate functions x M 
are the usual Einstein-Lorentz ones and have a precise operational meaning: x° = ct^\ 
where t is measured by "ideal clocks" at rest on I and synchronized "a la Einstein", 
x l , i = 1,2,3 are determined with ideal rules Jl3|]. (We use units where c = 1.) 

Let e M G sec TM be an orthonormal basis g(e fl ,e l/ ) = T]^ and e M = d/dx^^i, v = 
0,1,2,3). e determines an IRF. Let 7^ S sec/\ 2 (M) C secC€(M) be the dual basis 
and let 7 M = t] VlV 'f be the reciprocal basis to 7^, i.e., 7^.7^ = 6%. We have 7 M = da; M . 

As is well known the electromagnetic field is represented by a two-form F € sec /\ (M) C 
sec C£(M). We have 



/ -i; 1 -s 2 -i; 3 \ 

E 1 -B 3 B 2 

£; 2 B 3 -B 1 

\ £; 3 -B 2 B 1 / 



(2.72) 



where (E 1 , E 2 , E 3 ) and (B 1 ,^? 2 ,^ 3 ) are respectively the Cartesian components of the 
electric and magnetic fields. 

Let J e sec /^(M) C sec C^(M) be such that 



J = J^lp, = P7o + J 1 !! + J 2 12 + J 3 73, 



(2.73) 



where p and (J 1 , J 2 , J 3 ) are the Cartesian components of the cha rge a nd (3-dimen siona l) 
current densities. Recalling the definition of the operators d (eq.( 2.67 )) and 6 (eq.( 2.68| )) 
we see that we can write Maxwell equations as 



dF = 0, 5F= - J. 



(2.74) 



Since dF and SF are sections of C£(M) we can add the two equations in eq.( 2.74 ) 
and get 



(d - S)F = J. 



But from eq.( 2.69| ), d — S = d, the Dirac operator acting on sections of C£(M), and 
we get 

dF = J (2.75) 

which may now be called Maxwell equation, instead of Maxwell equations. 

We now write Maxwell equation in C£ + (M), the even sub-algebra of C£(M). The 
typical fiber of C£ + (M), which is a vector bundle, is isomorphic to the Pauli algebra 
(see Section 2.4). 



8 Here c is the constant called velocity of light in vacuum. In view of the superluminal and subluminal 
solutions of Maxwell equations found in this paper we don't think the terminology to be still satisfactory. 
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We put 

ffi = Jijo, i = = 7o7i7273 = 75- (2.76) 

Recall that i commutes with bi vecto rs and sinc e i 2 = — 1 it acts like the imaginary unit 
i = \^-T in C( + {M). From eq.( |2.69 ), using eq.( |2.35 ) we get 

F = E + iB (2.77) 

with E = E l 3i, B = B j aj, i,j = 1,2, 3. 

Now, since d = y^d" we get dy = d/dx° + 5^ = d/dx° - V. Multiplying eq.&7%) 
on the right by 70 we have 

<97o7o^7o = Jjo, 
(d/dx° — V)(— -B + i-B) = p + J, (2.78) 

where we us ed 7 °-F7o = —E + iB and J — J l c?;. 
From eq.( 2.78 ) we have 

-d E + id B +V£- iVB = p+J 

-d E + id B + X7.E + V A E - iV.B -iV AB = p + J 

Now we have 

-iVAiEVxi (2.79) 

since the usual vector product between two vectors a = a l c?i, b — b l di can be identified 
with the dual of the bivector a A b through the formula ax b = -i(aAi). Observe that 
in this formalism a x b is a true vector and not the nonsense pseudo vector of the Gibbs 
vector calculus. Using eq.(2.76) and equating the terms with the same grade we have 

V.E = p; VxB - d E = J ; 

(2.80) 

VxE + d B = 0; V.B = 0; 

which are Maxwell equations in the usual vector notation. 

We now introduce the concept of Hertz potential |iif| which permits us to find non- 
trivial solutions of the free "vacuum" Maxwell equation 

dF = (2.81) 

once we know nontrivial solutions of the scalar wave equation, 

□ $ = (d 2 /dt 2 - V 2 )$ = 0; $ e sec A°(M) c secCl(M) . (2.82) 

Let A E sec/\ 1 (M) c sec C£(M) be the vector potential. We fix the Lorentz gauge, 
i.e., d.A = -5A = such that F = dA = (d - 5) A = dA. We have the following 

Theorem. Let ir g sec/\ 2 (M) C secC£(M) be the so called Hertz potential. If ir 
satisfies the wave equation, i.e., On — d 2 n = (d — S)(d — 5)ir = — (dd + Sd)ir = and if 
we take A = —Sir, then F = dA satisfies the Maxwell equation dF = 0. 
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Proof. A = -Sir implies that 5A = -5 2 n = and F = dA = dA. Then dF = 
(d - 8)(d - 6) A = 5d(Sn) = -6 2 dir = 0, since Sdn = -dSi: from d 2 n = 0. 

From the above we see that if $ € sec f\°(M) c secC£{M) satisfies = 0, then 
we can find a non trivial solution of dF — 0, using a Hertz potential given, e.g., by 

7T = $7i72 ■ (2.83) 



In Section 4.3 this equation is used, e.g. to generate the superluminal electromagnetic 
JT-wave. 

We now express the Hertz potential and its relation with the E and B fields, in order 
for our reader to see more familiar formulas. We write it as sum of electric and magnetic 
parts, i.e., 

7T = 7r e + in m 



TT e = -^ % d; L , 7T m = -7T 23 CTl + 7T 13 <7 2 - 7r 12 CT 3 

Then, since A = dir we have 

A = ^(dir-nd) 

A'jQ = -3 7f e + V.7? e - (V X 7T m ) 

and since A = A^j^ we also have 

d_ 



(2.84) 



A = V.7? e ; A = A l cji = -t^tttc -Vxif„ 



d -» 

Since £ = -VA° - tt-^A B = Vx Awe obtain 
dx v 

E = -<9 (V x 7r TO ) + V x V x 7? e ; 

B = Vx (—do7T e - V X 7T TO ) = — <9o(V X 7T e ) - V X V X 7T m . 

We define E e ,B e , E m ,B m by 

£ e =VxVx7r e ; B e = -9 (Vxjr e ); 
£ m = -ob(V x 7r ro ) ; £?,„ = -V x V x 7? m . 

We now introduce the 1-forms of stress-energy. Since dF = we have Fd = 0. 
Multiplying the first equation on the left by F and the second on the right by F and 
summing we have: 

l/2(FdF + FdF) = d^{(l/2)F^F) = d^T" = 0, (2.86) 

where Fd = -{d^\F af}1 a ^)^ . 
Now, 

-J(W = "J(W) (2.87) 
Since j^.F = -(^F - FY) = F.^, we have 



2 



= -((F.<y*)F<f) i( 7 ^V)o 

= -(F.^).(FY) ~ \{F.Fyf.-f (2.88) 
= F^F x ^ aX + ^F a/3 F a P, 
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which we recognize as the stress-energy momentum tensor of the electromagnetic field, 
and = T» v ~f v . 

By writing F = E + \B as before we can immediately verify that 



T = -\f 1o F 



(E 2 + B 2 ) + (Ex B) 



7o- 



We have already shown that 9 M T M = 0, and we can easily show that 

d.T^ = 0. 



(2- 



(2.90) 



We now define the density of angular momentum. Choose as before a Lorentzian chart 
of the maximal atlas of M and consider the 1-form x = x^-f^ — x^^ . Define 

M„ = x A T M = i ix a T^ v - x,T a ^)j a A Y . 
It is trivial to verify that as = T vfl and d )1 T^ 1 ' = 0, it holds 



o. 



The invariants of the electromagnetic field F are F.F, F A F and F 
with 



F.F = --F^F^ ■ F A F = -J 5 F» V F^e^p 



(2.91) 
F.F + F AF 

(2.92) 



Writing as before F = E + \B we have 

F 2 = (E 2 - B 2 ) + 21E.B = F.F + F A F. (2.93) 
2.6. Dirac theory in C£{M) 

Let S = {7^} G sccA^M) C sec C£(M) be an orthonormal basis. Lett's € sec(/\°(M)+ 
/\ 2 (M) + /\ 4 (M)) C secC£(M) be the representative of a Dirac-Hestenes Spinor field 
in the basis E. Then the representative of Dirac equation in C£(M) is the following 
equation (H = c = 1): 

<9^£7i72 + rmps'-fo — . (2.94) 

To see that, consider the complexification C£c(M) of C£{M) calle d th e complex Clifford 
bundle. Then C£ c (M) = <U ® C£(M) and by the results of Section we know that the 
typical fiber of C£ C (M) is C4,i = (F ® C4,3, the Dirac algebra. 

Now let {ro,ri,r 2 ,r3,r 4 } C sec/\ 1 (M) C sec C£c(^) be an orthonormal basis with 



r a r& + r;,r a = 2^6 , 

9ab = diag(+l, +1, +1, +1, -1) . 



(2.95) 



Let us identify 7^ = T^T^ and call I = T0T1T2T3T4. Since I 2 — —1 and / commutes 
with all elements of C£^\ we identify / with i = \J — 1 and 7^ with the fundamental set 
of C£(M). Then if A G secC£ c (M) we have 



(2.96) 
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where $ s , $ p , A£, B£ u , t^ p € sec <U ® f\°(M) c secCl c (M), i.e., Vx € M, 
$ p (x), Ag(x), B^(x), Tg Vp {x) are complex numbers. 
Now, 

/= + 7o)J(l + i7i72) 



is a pr imiti ve idempotent field of C£c(M). We recall that, by eg. ( 2. 51 ), i/ = 7271/. 
From ( 2.91 ) we can write the following equation in C^c(-^0 : 

<9t/>s7i7 2 / + mi/>s7o/ = 
di/i^if — mipsf = 



and we have the following equation for \E' = "0s / : 

icM - = 0. 



(2.97) 



By eq.(2.49) and using for 7^ the matrix representation eq.(2.47) (denoted here by 
7 ) we get that the matrix representation of eq. ( 2.94| ) is 



ij^d^} -m|*) = 



(2.98) 



where now |$) is a usual Dirac spinor field. 

We now define a potential for the Dirac-Hestenes field tps ■ Since ips £ sec C£ + (M) it 
is clear that there exist A, B £ sec/\ (M) C sec Cl(M) such that 



=d(A + l5 B), 



(2.99) 



since 



d(A + 75S) = 9. A + d f\ A — lb d.B - j 5 d A B 
= S + B + l5 P- 

S = d.A; B = d l\ A — 75$ A B; P = -d.B . 
We see that when m — 0, satisfies the Weyl equation^ 

aVs = . 



Using eq.(2.10C) we see that 



d 2 A = d 2 B = 0. 



(2.100) 



(2.101) 



This last equation allows us to find UPWs solutions of arbitrary speeds for the Weyl 
equation once we know UPWs solutions of the scalar wave equation □ $ = 0, $ 6 
sec/\°(M) c secC£(M). Indeed it is enough to put A = (A + 75B) = $(1 + 75)1;, where 
v is a constant 1-form field. This result has been used in C^] to present subluminal 
and superluminal solutions of the Weyl equation. An example of a subluminal solution 
(indeed a stationary one) of the massless Dirac equation is obtained with the use of the 
"superpotential" Aq: 



C 

Ao = — (sin Sir cos Cit'-f — sinSlr sinSl,t'y 1 'y 2 j 3 ). 
r 



(2.102) 



9 We recall again that a Weyl spinor must satisfy 751/) = ±^721 (see e.g. ]10(]). 
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We have then 



C 

tpo = dAo = — [— fir 2 sin Sir sin Sli 

+ j° Xx cos fit + j°j 2 Xy cos fit 

+ 7°7 3 AzcosrM - 7 1 7 2 Azsinf2t 

+ 7 1 7 3 A?/sinrii - 7 2 7 3 Aa;sinfit 

+ 7°7 1 7 2 7 3 rir 2 sin fir cos fit], 



(2.103) 



where A = Sir cos Sir — sin Sir, r = \J x 2 + y 2 + z 2 . 
The above solution in the usual formalism reads 



. ' Xz .SI . 
lSinSn | + 1 — sin sir 

r>0 rv> 



*0 = 



i sin fit 
cos fit 



x 



\ 



r 

cos fit 



Xz 

~~3 



■ n • O 
1 — sin sir 

r 



A 



(2.104) 



/ 



Another very interesting possibility for constructing solutions of Weyl equation is the 
following. Suppose that F £ sec/\ 2 (M) C secC^(Af), F 2 ^ 0, is a solution of dF = 0. 
Then ip = e F is a massless Dirac spinor field satisfying dip = 0j^| Using this result and 
eq.(2.61) we can construct solutions of dipw = 0. To end this Section we show how to 
construct luminal or superluminal solutions of the Dirac equation. 

We know (see Section ||) that the Klein-Gordon equation has besides the subluminal 
solutions also luminal and superluminal solutions. Let $ be a subluminal, luminal or 
superluminal solution of + m$ = 0. Suppose $ is a section of C£c(M). Then in 
C£c{M) we have the following factorization: 



(d + \m)(d — im)<f> = 0. 

Now 

* = (d - im)$/ 
is a Dirac spinor field in C£c(M), since 

(9 + irn)* = 0. 



(2.105) 



(2.106) 



(2.107) 



'J is then a subluminal, luminal or superluminal UPW solution of Dirac equation, 
depending on <3>. 

3. EXTRAORDINARY SOLUTIONS OF THE (SCALAR) HOMOGENEOUS WAVE 
EQUATION AND OF KLEIN-GORDON EQUATION 



in 



A closed expression for e is given in ]16|. In particular, writing F = E + iB = zF (eq.(2.77)) 



where F = z*F/\z\ 2 , F 2 = 1, z 6 C, and F is a complex vector (in the Pauli algebra sense), then 
e F = e zF = cosh z + F sinh z. 
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3.1. Subluminal and superluminal solutions of the HWE 

Consider the HWE (c = 1) 

|!r$-V 2 4> = 0. (3.1) 

We now present some subluminal and superluminal solutions of this equation. 

Subluminal and Superluminal Spherical Bessel Beams. To introduce these beams we 
define the variables 

U^^ + y' + lKz-v^f]^; (3.2) 

(> = hx 2_ y 2 +1 2 >{z _ v>tr] l/2. (34) 

7> = — f== ; ^> - fc> = -0> ; w>=dcj > /dfc > . (3.5) 

v«> - 1 

We can now easily verify that the functions <&< m and below are respectively 
subluminal and superluminal solutions of the HWE (see example 3 below for how to 
obtain these solutions). We have 

$l m (t, x) = C e je(n p ^)P r e n (cos otf"* jl***-*,') (3.6) 

where the index p =<, >, Ci are constants, ji are the spherical Bessel functions, 
are the Legendre functions and (r, 6, ip) are the usual spherical coordinates. Q l ™ [^^J 1 ] 
has phase velocity (w < /k < ) < 1 [{w > /k > ) > 1] and the modulation function J^(^<£<) 
\je(fl>£>)] moves with group velocity u< [u>], where < v < < 1 [1 < u> < oo]. Both 
and <E>^ ra are undistorted progressive waves (UPWs). This term has been introduced 
by Courant and Hilbert Q|. However, they didn't suspect of UPWs moving with speeds 
greater than c = 1. For use in the main text we write the explicit form of <&<° and <I>2, , 
which we denote simply by <&< and $>: 

%(t,x) = C S[n( -^ p) e'K*-*") ; p=< or > . (3.7) 

When v< = 0, we have $< — > $o, 

£) = C £!!^ e iO< t; r = (a .2 + y 2 + z2) l/2 _ (3 8) 

r 

This solution has been found by Bateman in 1915 [0. The superluminal solution $> 
was discovered by Barut and Chandola in 1993 E^] . In what follows we show methods to 
obtain the solutions <!>< and $> for the HWE. When w> = oo, lo > = and $2, — > $00, 

* o(*,5)=a»^e <n >«, p^ + y 2 ) 1 / 2 . (3.9) 

We observe that if our interpretation of phase and group velocities is correct, then 
there must be a Lorentz frame where <fr p is at rest. It is trivial to verify that in 
the coordinate chart (x M ) which is a (nacs//') (see Section || for more details), where 
I' = (1 — 1> 2 ) -1 / 2 <9/9£+ (i></-\/l — v\)d/dz is a Lorentz frame moving with speed v < in 



the z direction relative to I = d/dt, $ p goes in $o(t' , x') given by eq.(3.8) with 1 1— > t' 
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x t— > x 1 . We can also verify that there is no Lorentz frame with velocity parameter 
< v < 1 where $> is at rest. 

Subluminal and Superluminal Bessel Beams. The solutions of the HWE in cylindrical 
coordinates are well known p[ . Here we recall how these solutions are obtained in order 
to present new subluminal and superluminal solutions of the HWE. In what follows the 
cylindrical coordinate functions are denoted by (p,9,z), p — (x 2 + y 2 ) 1 ^ 2 , x = pcos#, 
y = p sin 9. We write for $: 

3(t,p,M) = /i(p)/2(0)/s(t,*). (3-10) 



Inserting ( p.l0| ) in (|3.l| ) gives 



P 2 ^h + Pj-h + {Bp 2 - v 2 )h = (3.11) 

^+- 2 )/ 2 = 0; (3.12) 
d 2 d 2 \ 

^-^+^/ 3 = 0. (3.13) 

B and v are separation constants. Since we want $ to be periodic in 9 we choose v = n 
an integer. For B we consider two cases: (i) Subluminal Bessel solution, B = > 

In this case (3.11) is a Bessel equation and we have 

<S>%{t,p,9,z)=C n J n ( P n < y ( - k < z ~ w < t+ne \ n = 0,1,2,..., (3.14) 

where C n is a constant, J„ is the n-th order Bessel function and 

u}\ - k\ = n% . (3.15) 

In jl9) the <&j are called the n-th order non-diffracting Bessel beams.|^| Bessel beams 
are examples of undistorted progressive waves (UPWs). They are "subluminal" waves. 
Indeed, the group velocity for each wave is 

v< = dw</dfc<, < v< < 1 , (3.16) 

but the phase velocity of the wave is (aj</fc<) > 1. That this interpretation is correct 
follows from the results of the acoustic experiment described in J^, . 

It is convenient for what follows to define the variable rj, called the axicon angle [ po[ . 

fc< = k< cos 77 , f2< = k < sin 77 , < rj < tt/2 . (3-17) 

Then 

fc<=cj<>0 (3.18) 
„ = <J>^ 



and eq.(3.14) can be rewritten as <I>^ = <i>j , with 



$< = C„J„(/c<psin77)e i(fc < zcos,) ""<* + ™ e) . (3.19) 



11 The only difference is that k < is denoted by /3 = W — and w< is denoted by k' = u)/c > 0. 
(We use units where c = 1.) 
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9[ the n-th order non-diffracting portion of the 



In this form the solution is called in 
Axicon Beam. 

Now, the phase velocity v ph = 1/ cos 77 is independent of k 
dependent on k < 
can be subluminal or superluminal ! 

(ii) Superluminal (Modified) Bessel Solution, B = — f2> < 



but, of course, it is 

We shall show below that waves constructed from the 4>j beams 



In this case (3.11) is the modified Bessel equation and we denote the solutions by 
(t,p,9,z) =CA(f!>p)e i(i >" u >t+n9 \ n = 0,l,..., (3.20) 



where K n are the modified Bessel functions, C n are constants and 



(3.21) 



We see that are also examples of UPWs, each of which has group velocity w> = 
dw > /dfc> such that 1 < v > < 00 and phase velocity (oj>/fc>) < 1. As in the case of 
the spherical Bessel beam (eq.(3.7)) we see again that our interpretation of phase and 
group velocities is correct. Indeed, for the superluminal (modified) Bessel beam there 
is no Lorentz frame where the wave is stationary. 

The &x beam was discussed by Band |2^] in 1988 as an example of superluminal 
motion. Band proposed to launch the <&^- o beam in the exterior of a cylinder of radius 
r\ on which there is an appropriate superficial charge density. Since iio(^> r i) is non 
singular, his solution works. In Section ^ we discuss some of Band's statements. 

We are now prepared to present some other very interesting solutions of the HWE, 
in particular the so called X-waves, which are superluminal, as proved by the acoustical 
experiments described in |(| . 

Theorem (Lu and Greenleaf ) The three functions below are families of exact solu- 
tions of the HWE [eq.(3.1 )] in cylindrical coordinates: 



= / T(fc<) 



^ K ( S ) = / D(r,) 



1 

2T 
1 

2^ 



A(4>)f(a)d4> 



A(4>)f{s)c\4> 



dfc< 
dr) ; 



where 
and 



$ L (p,6,z-t) = $i(/?,0)$ 2 (2-t) ; 
s = ao(k < ,f])pcoa(8 — 0) + b(k < ,n)[z ± Ci(fc<, rj)t] 



ci(k<,v) = V 1 + [ao(k <1 ri)/b(k <1 r))} 2 



(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 



In these formulas T(k < ) is any complex function (well behaved) of k< and could 
include the temporal frequency transfer function of a radiator system, A(cf>) is any com- 
plex function (well behaved) of <f> and represents a weighting function of the integration 
with respect to 4>, f(s) is any complex function (well behaved) of s (solution of eq.(3.1)), 
D(rj) is any complex function (well behaved) of 77 and represents a weighting function 
of the integration with respect to 77, called the axicon angle, ao(A<,))) is any complex 
function of k < and rj, b{k < ,rf) is any complex function of k < and r\. 
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As in the previous solutions, we take c = 1. Note that k < , rj and the wave vector k < 
of the f(s) solution of eq.(|3.lD are related by eg. ( [3. 17 ). Also $2(2 — t) is any complex 
function of (z — t) and ^\{p, 9) is any solution of the transverse Laplace equation, i.e., 



"1 d 




1 d 2 ' 


pdp 




+ p^W 



3i(p,0) = O. 



(3.27) 



The proof is obtained by direct substitution of (f> r) , $^ and in the HWE. Obviously, 
the exact solution is an example of a luminal UPW, because if one "travels" with 
the speed c = 1, i.e., with z — t = constant, both the lateral and axial components, 
<l>i( p, 8) and $2(z — t) will be the same for all time t and distance z. When ci(k,rf) in 
eq.( 3.25| ) is real, (±) represent respectively backward and forward propagating waves. 

We recall that & v (s) and $k(s) represent families of UPWs if ci (&<,??) is indepen- 
dent of k < and 77 respectively. These waves travel to infinity at speed c%. $jj(s) is a 
generalized function that contains some of the UPWs solutions of the HWE derived 
previously. In particular, if T(/«<) = S(k < — fc < ), k < = uj > is a constant and if 
/(s ) = c s , ao(k < ,rj) = — if2<, b(k < ,rf) = \j3 = iw/ci, one obtains Durnin's UPW beam 



Durnin 



-if2<pcos(# — 0) 



iiPx-ut) 



(3.28) 



-th order UPW Bessel beam $j given by eq.(3.14) 



If A(<£) = i n e m *, we obtain the n- 
(s) is obtained in the same way with the transformation k < = k < cosrj; fl < = 
k < sin 77. 

The X-waves. We now present a superluminal UPW discovered in 1992 by Lu and 
Greenleaf [[[9) which, as discussed in [|[ ^(|, is physically realizable in an approximate 
way in the acoustic case and can be used to generate Hertz potentials for the electro- 
magnetic field (see Section 0). We take in eq.(3.22): 



T(fe<) = 5(fc<)e~ a ° fc < ; A(<f>) = i n e in *; /(«) = e s 
ao(fc<, 77) = — ifc< sin 77 b{k < , 77) — ifc< cos 77; . 



(3.29) 



Then we get 



$ 1 



B(fc< ) J„ (jfc<p sin ?7)e- fc < [ao cos "-* )] dfc< 



(3.30) 



In eq.(3.30) B(k < ) is any well behaved complex function of fc< and represents a transfer 
function of a practical ra diato r, k < — lj and ao is a constant, and 77 is again called the 
axicon angle. Equation ( |3.30|) shows that ^> x is represented by a Laplace transform 
of the function B(k < ) J n (k < p sin 77), and an azimuthal phase term e . The name X- 
waves for the $5^ comes from the fact that these waves have an X-like shape in a plane 
containing the axis of symmetry of the waves ||, ||] . 



The ^xbb^ waves. This wave is obtained from eq.( |3.30| ) putting B(k < ) = do. It is 
called the X-wave produced by an infinite aperture and broad bandwidth. We use in 
this case the notation & X bb • Under these conditions we get 



ao (p sin77) n e 



<M{t + VM) 



(n= 0,1,2,...) 



(3.31) 
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where the subscript denotes "broadband" . Also 

M = (psuiTy) 2 +t 2 ; (3.32) 
r = [do — i(z cos 77 — t)] (3.33) 

For n = we get $j BB and 

" Ai ' Jj " ^(psin^) 2 + [a - i(zcos7/ - t)] 2 



It is clear that all $^ BB are UPWs which propagate with speed c\ = 1/ cosry > 1 in 
the ^-direction. Our statement is justified for as can be easily seen (as in the modified 
super luminal Bessel beam) there is no Lorentz frame where &xbb ^ s a * res t- Observe 
that this is the real speed of th e wave; phase and group velocity concepts are not 
applicable here. Equation ( [3.34 ) does not give any dispersion relation. Th e $y RR 
waves cannot be produced in practice as they have infinite energy (see Section 3.4), but 
a good approximation for them has been realized with finite aperture radiators |)23|, ^| . 

Recall that if in cq.( 3.30 ) we put B(k<)e~ a ° k< — A(k < ) and if we take into account 
that for each component Bessel beam in the packet the following dispersion relation 
holds: 

lu = k < / cosrj = k/ cosr] (fc< = k) , (3.35) 

then for both the broad band X-waves (as e.g. eq.( |3.34[ )) and the limited band X-waves 
where A(k) = A{k < ) is centered in ko we can write 



w{k) = uj(k ) + ^ 



(k-k ) = Lu + ^(k-k ), (3.36) 



where coq — 1/ cosr/, dui/dk — 1/ cos rj. It follows that eq.(3.30) can be written 

$x (t,x,y,z) = / dkA(k) JJkpta,nr])e i{z ~~ t)k 

cos 7] J 

= $(o, x ,y, x -—t) (3.37) 
\ cos 77 / 

showing that the X-waves propagate without distortion with speed dw/dko = 1/ cost?. 

We end this section with the commentary that in j^] we develop methods for project- 
ing "finite aperture approximations" to the exact acoustic and electromagnetic solutions 
discussed in this paper. 

3.2. Donnelly- Ziolkowski method for designing subluminal, luminal and su- 
perluminal UPW solutions of the HWE and the Klein-Gordon equation 
(KGE) @ [25 1 

Consider first the HWE for $ (eq.( |3.l| )) in a homogeneous medium. Let $(w, k) be the 
Fourier transform of <!>(£, x), i.e., 

$(w,fc)= / d 3 x dtm,^-'^-^, (3.38) 

JR 3 J-oo 

1 f -> r + °° ~ 
$(t,x) = -—r d 3 k duj^ioj,^ 1 ^-^. (3.39) 



(2tt)' 
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Inserting J3.38D in the HWE we get 

{uj 2 - k 2 )$>{uj,k) = (3.40) 

and we are going to look for solutions of the HWE and eq.( 3.40| ) in the sense of distri- 
butions. We rewrite eq.(3.40) as 

{uj 2 - k 2 z - 2 )$(w,fc) = 0. (3.41) 

It is then obvious that any $(w, k) of the form 

$(u>, k) = 5(0, /3) 5[uj - (/? + 2 /4/?)] 6[k z - {p - 2 /4/?)] , (3.42) 

where 5(0,/?) is an arbitrary weighting function, is a solution of cq.(3.41) since the 
5-functions imply that 

uj 2 - k\ = Q 2 . (3.43) 

In 1985 Ziolkowski |2(J found a luminal solution of the HWE called the Focus Wave 
Mode. To obtain this solution we choose, e.g., 

E FWM {£l,P) = % exp(-O 2 z /4/3), (3.44) 

whence we get, assuming (3 > and zq > 0, 

9 FWM {t,x) - e yy — — — . (3.45) 

Despite the velocities v% = +1 and V2 = —1 appearing in the phase, the modulation 
function of &fwm has ver y inte resting properties, as discussed in details in ^6|. It 
remains to observe that eq.( |3.45 ) is a special case of Brittingham's formula p7[ . 

Returning to eq.( 3.42| ) we see that the (5-functions make any function of the Fourier 
transform variables uj, k z and fi to lie in a line on the surface w 2 — k 2 — Vt 2 = (eq.( |3.41 )). 
Then, the support of the (5-functions is the line 

lu = f3 + n 2 /4(3; k z = /3- n 2 /A(3 . (3.46) 

The projection of this line in the (to, k z ) plane is a straight line of slope —1 ending at 
the point (/?, j3). When (3 = we must have = 0, and in this case the line is uj = k z 
and <I>(i, x) is simply a superposition of plane waves, each one having frequency uj and 
traveling with speed c = 1 in the positive z direction. 

Luminal UPWs solutions can be easily constructed by the ZM, but will not be dis- 
cussed here. Instead, we now show how to use ZM to construct subluminal and super- 
luminal solutions of the HWE. 

First Example: Reconstruction of the subluminal Bessel Beams $ j and the superlu- 
minal &xbb (^-wa ve )- Starting from the "dispersion relation" uj 2 — k 2 z — fl 2 = 0, we 
define 

$(w,fe) = E(k,ri)5(k z - k cos r])S (uj ~k). (3.47) 

This implies that 

k z = kcosrj; cosrj = k z /uj, uj > 0, — l<cos?y<l. (3.48) 
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Wc take moreover 



£1 = A; sin 77; k > 



(3.49) 



We re call th at ft — (k x ,k y ), p — (x,y) and we choose VI. p = VLpcosO. Now, putting 
eq.( 3.47 ) in eq.( 3.39| ) we get 



x) 



Choosing 



(2-) 4 J 



dk fcsin rj 



d6E(k,rj)e l 



kp sin rj cos 6 



i(fc cos rjz — kt) 



~(M) = (2.) 3 ^ 



kzo sin rj 



k sin 77 



where zq > is a constant, we obtain 



3>(i, a;) = zq sin 77 



1 

2^ 



2tt 



d$ e^^ s * n ^ cos ^ 



Calling zosin?7 = a-o > 0, the last equation becomes 

roc _ _ 

$> o (M) = a / dIe" feao J (fcpsin?7)e i ' c(cos? ' 2 - t) 



ifc(cos 7] z — t) 



Writing fc = fc< and taking into account eq.(3.19) we see that 

J (fc<psin77)e iI <^ cos "-*) 



(3.50) 



(3.51) 



(3.52) 



(3.53) 



(3.54) 



is a subluminal Bessel beam, a solution of the HWE moving in the positive z direction. 
Moreover, a comparison of eq.( 3.53 ) with eq.( 3.30 ) shows that ( 3.53 ) is a particular 
sup erlum inal X-wave, with B(k < ) — e~ a ° k< . In fact it is the $xbb UPW given by 
eq.(l3~34l). 

Second Example: Choosing in (3.50) 

H(fc,r?) = (27r) 3 e- 2o|cosii|I cot7/ (3.55) 

gives 



<S>>(t,x) 



ros'r?/ dfc/fce- 2o|cosr ' l ' £ Jo(fcpsinr/)e~ i ' £(cos " z ~ t) (3.56) 
[20 — isgn(cosry)(z — tj cost?)] 



[p 2 tan 2 rj + [zq + isgn(cos7])(z — tj cos?7)] 2 ] 3 / 2 



(3.57) 



Comparing eq.(3.56) with eq.( |3.30 ) we discover that ZM produced in this example 



a more general $^ wave where -B(fc<) 



-zq\ cos rj\k < 



Obviously $ > (t,a;) given by 



cq.(3.57) moves with superluminal speed (l/cosr;) in the positive or negative z-direction 
depending on the sign of cos 77, denoted sgn(cosry). 

In both examples studied above we see that the projection of the supporting line of 
eq.(3.47) in the (w, k z ) plane is the straight line k z /oj = cos 77, and cos 77 is its reciprocal 
slope. This line is inside the "light cone" in the {uj, k z ) plane. 
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Third Example: Consider two arbitrary lines with the same reciprocal slope that we 
denote by v > 1, both running between the lines u = ±k z in the upper half plane uj > 
and each cutting the w-axis at different values (3i and /3 2 . 

The two lines are projections of members of a family of HWE solution lines and each 
one can be represented as a portion of the straight lines (between the lines to = ±fc 2 ) 



v(u> - /3i), k z = v(u - /3 2 ). 



(3.58) 



It is clear that on the solution line of the HWE, £1 takes values from zero up to a 
maximum value that depends on v and (3 and then back to zero. 

We see also that the maximum value of il, given by j3v/y/v 2 — 1, on any HWE 
solution line occurs for those values of to and k z where the corresponding projection 
lines cut the line u = vk z . It is clear that there are two points on any HWE solution 
line with the same value of il in the interval 



< Q < vf3/y/v 2 -1 = Q . 



(3.59) 



It follows that in this case the HWE solution line breaks into two segments, as is the 
case of the projection lines. We can then associate two different weighting functions, 
one for each segment. We write 



$(n,w,k z ) = 3i(fi,«,/3)5 



vip + ^/^v 2 - n 2 (v 2 - 1)] 

(v 2 - 1) 



x5 



[fiv 2 + ^Jv 2 l3 2 -Q 2 (v 2 -1)] 



(v 2 - 1) 

V[(3- y/p 2 V 2 -fl 2 (v 2 -1)} 

(v 2 - 1) 



[f3v 2 - ^Jv 2 f3 2 - Q 2 (v 2 - 1)] 
(v 2 - 1) 



Now, choosing 



we get 

$„ i/3 (£,p, z) = f2 exp 



i/3v(z — vt) * 



n 2 



(3.60) 



(3.61) 



dX XJ (Q, pX) cos 



ilov (z — t/v) 

Vv 2 ~ 1 y/l - X 2 



Then 



$ Vt p(t,p,z) = exp 



v(z — vt) 

7^ 



1 



in {n ^j^(z-t/v) 2 +p 2 } 
{^^{z-t/vY+p 2 ) 



(3.62) 
(3.63) 



If we call v < — — < 1 and take into account the value of Qq given by eq.(3.59), we can 

r- — . v ' v 

write eq.( 3.63 ) as 



<S> v< (t,P,z) = S[ni ^ <] e in ^-^ ; 



(3.64) 
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x 2 + y 2 + 



1 



1 - V' 



-{z- V< tf 



1/2 



(3.65) 



which we recognize as the subluminal spherical Bessel beam of Section 3.1 (eq.(p.7|)). 
3.3. Klein-Gordon equation (KGE) 

We show here the existence of subluminal, luminal and supcrluminal UPW solutions of 
the KGE. We want to solve 



d 2 

_ - V 2 + m" ) $ K<J (t,f) = 0, m>0, 



(3.66) 



with the Fourier transform method. We obtain for $ kg (uj, k) (a generalized function) 
the equation 

{uj 2 -k 2 z - (n 2 + to 2 )}<F g (w, k) = 0. (3.67) 

As in the case of the HWE, any solution of the KGE will have a transform $(w, k) 
such that its support line lies on the surface 



lo 2 - k 2 z - (Vl 2 + m 2 ) = 



(3.68) 



From eq.(3.68), calling ft 2 + m 2 = K 2 , we see that we are in a situation identical to 
the HWE for which we showed the existence of subluminal, supcrluminal and luminal 
solutions. We write down as examples one solution of each kind. 

Subluminal UPW solution of the KGE. To obtain this solution it is enough to change 



in eq.(3.63) n = vf3/Vv 2 - 1 -> fl, 



KG 




v/3 



1 1/2 



$£ G (i,p,z) =exp 



i/3v(z-vt) \ sin(f^ G £<) 



We have, 



x 2 + y 2 + 



1 



V^ - 1 



{z- V< tf 



-,1/2 



, 1>< = 1/v. 



(3.69) 



(3.70) 



Luminal UPW solution of the KGE. To obtain a solution of this type it is enough, 

(3.71) 



as in eq.(3.42), to write 

$ KG = S(fl, (3)S[kz - (n 2 + (to 2 - I3 2 )/2I3)]5[u - (ft 2 + (m 2 + (3 2 )/2(3)} 
Choosing 



gives 



B(n,j9) = ^cxp(-z r! 2 /2/3), z > 0, 



<f> K G(+ ^ _ „i Z (m 2 -/3 2 )/2/3^it(m 2 +/3)/2/3 ex P{"P 2 / 3 / 2 ^0 ~ K z ~ 

[z Q - l(z - t)\ 



(3.72) 



(3.73) 
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Superluminal UPW solution of the KGE. To obtain a solution of this kind we intro- 
duce a parameter v such that < v < 1 and write for <S> KG in ( |3.67 ) 



$Zjt(w,il,k.) = E(n,v,(3)5 



-(3v 2 + ^{Q 2 +m 2 )(l-v 2 )+v 2 /3 2 
1 - v 2 



x 8 



(-0 + v / (n 2 +m 2 )(l-v 2 )+v 2 P 2 ^j 
1-v 2 



Next we choose 



E(Q,v,0) 



(27r) 3 exp(-zoV fi o + ^ 2 ) 



where zo > is an arbitrary parameter, and where 

P 2 v 2 



i 



(3.74) 



(3.75) 



(3.76) 



Then introducing u> = 1/v > 1 and 7> 



1 



we get 



it^-m^U-rt) exp |-^o - i7>(^ - w>^)] 2 + + y 2 | 



(3.77) 



y/[z - i7>(z - Vyt)] 2 + x 2 + y 2 
which is a superluminal UPW solution of the KGE moving with speed v > in the z 



direction. From eq.(3.77) it is an easy task to reproduce the superluminal spherical 
Bessel beam which is solution of the HWE. 

3.4. On the energy of the UPWs and the velocity of transport of energy 

Let $ r (t,x) be a real solution of the HWE. Then, as it is well known ^8j, the energy 
of the solution is given by 



dv 



d<5> r 



$ r V 2 <I\ 



lim 



d5$ r n.V$ r 



(3.78) 



S(R) 



where S(R) is the 2-sphere of radius R . 

We can easily verify that the real or imaginary parts of all UPWs solutions of the 
HWE presented above have infinite energy. The question arises of how to project 
superluminal waves, solutions of the HWE, with finite energy. This can be done if we 
recall that all UPWs discussed above can be indexed by at least one parameter that 
here we call a. Then, calling $ a (t,x) the real or imaginary parts of a given UPW 
solution we may form "packets" of these solutions as 



J daF{a)<f> a (t,x) 



(3.79) 



We now may test for a given solution $ Q and for weighting function F(a) if the 
integral in eq.( 3.78| ) is convergent. We can explicitly show that for some (but not all) of 
the solutions showed above (subluminal, luminal and superluminal) that for weighting 
functions satisfying certain integrability conditions the energy e results finite. 
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It is particularly important in this context to quote that the finite aperture approxi- 
mations for all UPWs discussed in this paper have, of course, finite energy. For the case 
in which <I> given by eq.( 3.79| ) is used to generate solutions for, e.g., Maxwell of Dirac 
fields, the conditions for the energy of these fields to be finite will in general be different 
from the condition that gives for $ a finite energy. This problem will be discussed with 
more details in another paper. 

To finish we remark that for a scalar field satisfying 



1 92 n2 

V 

a dt 2 



$ = 



we have as is well known that the flux of momentum is given by 

5$ 



S = V$ 



and 



(V$) 2 



dt 

1 

cl \~dt 



(3.80) 



(3.81) 



(3.82) 



We can immediately verify that if the speed of transport of energy is defined as 
v E = I S\/u then v e < c*. The acoustic experiments reported in j^, |(| show nevertheless 
that for the X-wave the energy travels with speed c s / cos 77 (c* = c s ). We thus see that 
the usual definitions of magnitudes such as density of energy and momentum and the 
velocity of transport of e nergy demand a careful revision. (See in this context also the 
discussion of Section 4.4.) 



4. SUBLUMINAL AND SUPERLUMINAL UPW SOLUTIONS OF MAXWELL EQUA- 
TIONS (ME) 

In this Section we make full use of the Clifford bundle formalism (CBF) summarized in 
Section @, but we translate all the main results into the standard vector formalism used 



by physicists. We start by reanalyzing in Section 4.1 the plane wave solutions (PWS) 
of ME with the CBF. We clarify some misconceptions and explain the fundamental 
role of the duality operator 75 and the meaning of i = in standard formulations 

of electromagnetic theory. Next, in Section 4.2 we discuss subluminal UPWs solutions 
of ME and an unexpected relation between these solutions and the possible existence 
of purely electromagnetic particles (PEPs) envisaged by Einstein p9| , Poincare |50|, 
Ehrenfest |31 and recently discussed by Waite, Barut and Zeni ]32|, |3^|. In Section 
4.3 we discuss the theory of supcrluminal electromagnetic X-waves (SEXWs). In §, 
6| we present simulations of the motions of the SEXWs and of their finite aperture 
approximations, which can eventually be launched by appropriate physical devices. 



4.1. Plane wave solutions of Maxwell equations 



We recall that ME in vacuum can be written as [eq.(2.75)] 



OF = 0, 



(4.1) 



where F sec f\ (M) C sec Ci(M). The well known PWS of eq.([4.l|) are obtained as 
follows. We write in a given Lorentzian chart (x^) of the maximal atlas of M a PWS 
moving in the z-direction 

F = fe^ kx (4.2) 
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k = fc^, fc 1 = k 2 = 0, x = x^, (4.3) 

wher e fc, x € sec/\ (M) C secC£(M) and where / is a constant 2-form. From eqs.(4.1) 
and (4.2) we obtain 

kF = (4.4) 

Multiplying eq.(4.4) by fc we get 

k 2 F = (4.5) 
and since k G sec/\ 1 (M) C sec C£(M) then 

fc 2 = <-> fe = ±|fc| = k 3 , 
i.e., the propagation vector is light-like. Also 

F 2 = F.F + F A F = 



(4.6) 



(4.7) 



as can be easily seen by multiplying both members of eq.(4.4) by F and taking into 
account that k ^ 0. Eq(4.7) says that the field invariants are null. 

It is interesting to understand the fundamental role of the volume element 75 (duality 
operator) in electromagnetic theory. In particular since e 7sfc:E = coskx + 75sinfca;, we 
see that 

F = f cos kx + 75/ sin kx. (4-8) 

Writing F = E + iB, (see eq. fl2.72j )) with i = 75 and choosing / = e% + ie2, ei.e*2 = 0, 

becomes 



ei, e*2 constant vectors in the Pauli subalgebra sense, eq.([4 

E + iB = e*i cos kx — €2 sin kx + i(e*i sin kx + e*2 cos kx) 



(4.9) 



This equation is important because it shows that we must take care with the i = 
y/— 1 that appears in usual formulations of Maxwell theory using complex electric and 
magnetic fields. The i = \/— 1 in many cases unfolds a secret that can only be known 
through eq.(4.9). From eq.(4.4) we can also easily show that k.E = k.B — 0, i.e., PWS 
of ME are transverse waves. 
We can rewrite eq.(4.4) as 

fc7o7o^7o = (4.10) 
and since kjo = k$ + k, 70-F70 = — E + iB we have 

kf = k f. (4.11) 

Now, we recall that in Cl + (M) (where, as we said in Section 0, the typical fiber is 
isomorphic to the Pauli algebra €£3,0) we can introduce [Q the operator of space 
conjugation denoted by * such that writing / = e + ib we have 



f* = -e + ib 



ko = k 



k* = -k. 



We can now interpret the two solutions of k 2 — 0, i.e., fco — 1^1 & n d fco 



(4.12) 

— \k\ as 

corresponding to the solutions k f — kf and k f* = —kf*; f and /* correspond 
in quantum theory to "photons" which are of positive or negative helicities. We can 
interpret fco = 1^1 a s a particle and fco = — |fc| as an antiparticle. 

Summarizing we have the following important facts concerning PWS of ME: (i) the 
propagation vector is light-like, fc 2 = 0; (ii) the field invariants are null, F 2 — 0; (iii) 
the PWS are transverse waves, i.e., k.E = k.B = 0. 
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4.2. Subluminal solutions of Maxwell equations and purely electromagnetic 
particles 

We take $ S sec(/\°(M) © /\ 4 (Af)) C secC£(M) and consider the following Hertz 
potential it 6 sec A 2 (Af) C sec«(M) [eq.(p^)] 

7T = $7 X 7 2 - (4.13) 

We now write 

$(*,£) = 0(f)e 75f2t . (4.14) 
Since 7r satisfies the wave equation, we have 

V 2 0(f) + f2 2 0(x) = (4.15) 

Solutions of eq.(3.15) (the Helmholtz equation) are well known. Here we consider the 
simplest solution in spherical coordinates, 



r 

where C is an arbitrary real constant. From the results of Section || we obtain the 
following stationary electromagnetic field, which is at rest in the reference frame Z 
where are naturally adapted coordinates to Z (see Section || for the definition of 
these concepts). 

C 

F = — [sin Slt(a£lr sin 9 sin ip — (3 cos 9 sin 9 cos tp)^oli 

— sin Qt(aflr sin 9 cos tp + /3 sin cos sin y)7o72 

+ sin fii(/? sin 2 6> - 2a)7o7 3 + cos Vtt{(3 sin 2 - 2a)7i7 2 

+ cos fM(/3 sin cos sin (/? + afir sin cos <p)7i73 

+ cosftt(— /3 sin 6* cos # cos <^ + aOr sin sin ^3)7273] (4-17) 

with a = Or cosilr — sin Sir and (3 = 3a + il 2 r 2 sin fir. Observe that Fq is regular at 
the origin and vanishes at infinity. Let us rewrite the solution using the Pauli-algebra 
in C£+(M). Writing (i = 75) 

F =E + LB (4.18) 

we get 

E = WsinQt , B = WcoaQt (4.19) 

with 

/ afly 0xz aQx (3yz (3(x 2 + y 2 ) 2a" 



Or 



W = -C[-f-Z- r , 3--^, v 5 — ~— • (4-20) 

We verify that divW> = 0, div£ = divB = 0, mtE +dB /dt = 0, rotB -dE /dt = 0, 
and 

rotw = nw. (4.21) 
1 



Now, from eq.( 2.89 ) we know that To = — -F70-F is the 1-form representing the energy 



density and the Poynting vector. It follows that Eq x Bq = 0, i.e., the solution has zero 
angular momentum. The energy density u = S 00 is given by 

u = \ [sin 2 9{9?r 2 a 2 + cos 2 90 1 ) + {(3 sin 2 9 ~ 2a) 2 ]. (4.22) 
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Then J J j R3 udv — oo. As for the case of the scalar field (see Section 4.4) a finite 



energy solution can be constructed by considering "wave packets" with a distribution 
of intrinsic frequencies F(Q,) satisfying appropriate conditions. Many possibilities exist, 
but they will not be discussed here. Instead, we prefer to direct our attention to 
eq.(4.21). As it is well known, this is a very important equation (called the force free 
equation p^|) that appears e.g. in hydrodynamics and in several different situations in 
plasma physics pq]. The following considerations are more important. 



Einstein [|29| among others (see 1 32 for a review) studied the possibility of construct- 
ing PEPs. He started from Maxwell equations for a PEP configuration described by an 
electromagnetic field F p and a current density J p , where 



and rightly concluded that the condition for existence of PEPs is 



Jp.Fp — 0. 



This condition implies in vector notation 



PpE p = , jp.Ep = , j p x Bp = 



(4.23) 



(4.24) 



(4.25) 



From eq.(4.25) Einstein concluded that the only possible solution of eq.(4.23) with the 
subsidiary condition given by eq.(4.24) is J p — 0. However, this conclusion is correct, as 
pointed in 32, 3^|, only if > 0, i.e., if J p is a time-like current density. If we suppose 
that J p can be spacelike, i.e., J^ < 0, there exists a reference frame where p p = and 
a possible solution of eq.(3.24) is 



p p = , Ep.Bp = , j p = KCBp, (4.26) 
where K = ±1 is called the chirality of the solution and C is a real constant. In p2l p3 



static solutions of eq.( 4.23 ) and ( 4.24 ) are exhibited where E p = 0. In this case we can 
verify that B p satisfies 



VxB p = KCBp. 
Now, if we choose F £ sec/\ 2 (M) c secd(M) such that 

Fq = Eq + iB , 
Eq = Bp cos ilt , Bo = Bp sin fit 

and Q — KC > 0, we immediately realize that 

dF = . 



(4.27) 



(4.28) 



(4.29) 



This is an amazing result, since it means that the free Maxwell equations may have 
stationary solutions that model PEPs. In such solutions the structure of the field Fq is 
such that we can write 



F = F p + F = \W cos nt-W sin M , 



dF r , 



-OF = J, n 



(4.30) 



i.e., dFo — is equivalent to a field plus a current. This opens several interesting 
possibilities for modeling PEPs (see also ]^6)) and we discuss more this issue in another 
publication. 
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We observe that moving subluminal solutions of ME can be easily obtained choosing 
as Hertz potential, e.g., 



n < (t,x)=C 



sin f2£< 



exp[7 5 (w<i - k < z)]jij 2 , 



U = [x 2 + y 2 + 7 2(z - V< t) 2 ] , 



(4.31) 
(4.32) 



1 



7< 



v < = dw</dfe< 



We are not going to write explicitly the expression for F K corresponding to 7r < because 
it is very long and will not be used in what follows. 

We end this Section with the following observations: (i) In general for subluminal 
solutions of ME (SSME) the propagation vector satisfies an equation like eq.(4.30). (ii) 
As can be easily verified, for a SSME the field invariants are non-null, (iii) A SSME is 



not a transverse wave. This can be seen explicitly from eq.(4.20). 

Conditions (i), (ii), (iii) are in contrast with the case of the PWS of ME. In |37], Q 
Rodrigues and Vaz showed that for free electromagnetic fields (dF = 0) such that 
F 2 7^ 0, there exists a Dirac-Hestenes equation for ip £ sec(/\°(M)+/\ 2 (M)+/\ 4 (M)) C 
sec C£(M) where F — i/ , 7i72V'- This was the reason why Rodrigues and Vaz discovered 
subluminal and supcrluminal solutions of Maxwell equations (and also of Weyl equation 
p5[) which solve the Dirac-Hestenes equation [eq.(2.94)]. An explicit superluminal 
solution of Maxwell equations is given in jy| using as Hertz potential n> = $>7 12 
where <J>> is given by eq.(3.6). 



4.3. The superluminal electromagnetic A-wave (SEXW) 

In this Section we present a family of solutions of Maxwell equations called the super- 
luminal electromagnetic A-waves Fx BB n ■ A solution dual to Fx bb„ , called *Fx b B n 
has been first presented by Lu and Greenleaf in an unpublished paper ||| . Later the 
solutions Fx b b„ and others associated with it have been studied in detail [| . 

T o sim plify the matter in w hat follows we now suppose that the functions <I>x„ 
[cq.( 3.30| )1 and &xbb„ [eq.( 3.31 )1 which are superluminal solutions of the scalar wave 
equation are 0-fo rms sections of th e com plexi fied C lifford bundle C£c(M) = <S ' ®Cl(M) 
(see Section ^2~6|) . We rewrite eqs.( 3.30 ) and ( 3.34 ) as (n = 0, 1, 2, . . .) 



®x n (t, x) 



and choosing B(k) — qq we have 



B(fc)J n (fcpsin?7)e- fe[ao - i(2COST '- t) ]dA; 



<I>XBB n {t,x) = 
M = (psini]) 2 + t 2 



ao(p sin77) n e 1 ' 



<M(t + VM) n 

; r = [ao — i(z cos T) — t)]. 



(4.33) 

(4.34) 
(4.35) 



Further, we suppose now that the Hertz potential it, the vector potential A and the 
corresponding electromagnetic field F are appropriate sections of Cic(M). We take 



tt = $7172 G sec d7 ® f\ 2 (M) C secC£ c {M), 



(4.36) 
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where <£> can be &x„ , &xbb„ , &xbl„ ■ Let us start by giving the explicit form of the 
F X bb„ the SEXWs. In this case eq.(2.84) gives tt — i 7r m and 



= $xbb„z 



(4.37) 



where z is the versor of the z-axis. Also, let p, 6 be respectively the versors of the p 
and 9 directions where (p, 6, z) are the usual cylindrical coordinates. Writing 



FxBB n = ExBB n 
we obtain from equations ( 3.31 ) and ( [2.85 ): 

p d 2 



75-B 



XBB n 



E 



XBB n — 



pdtdO 



XBB n 



■6*-* 

dtdp 



XBB n 



D 



o 2 



XBB„ 



dpdz 



XBB„ 



e 1 - 



d 2 



p dOdz 



XBB„ 



dz 2 df 2 1 



Explicitly we get for the components in cylindrical coordinates, 
-» 1 M 

(E X BB n )p = --n—=®xBB„ 
P VM 

1. Me 

: -1 7=77 ®XBB n 

p \JMM% 



(ExBB n )e 

(BxBB n )p 

(BxBB n )e 

{BxBB n )z 



cosr](E X BB n )e 

- cosrj(E X BB n )p 
2 M 7 



sin i]- 



XBB„ 



The functions Mi(i = 2, 



IM 

, 7) in (4.41) are 



M 2 = T^ 
M 3 = n- 

M 4 = 2n 
M 5 = t H 



1 



M 
3 



M 6 
M 7 



2 „;„2 



(/? sin ?7 



(n 2 



1) 



M 
M 
M 4 

M~ 
1 



nM 3 )M 2 + np 



3n—T 

M 



2 M 5 ■ 2 

sm ii 

M 



3- 



1 



(4.38) 

(4.39) 
(4.40) 

(4.41a) 

(4.41b) 

(4.41c) 
(4.41d) 

(4.41e) 

(4.42a) 
(4.42b) 

(4.42c) 
(4.42d) 
(4.42e) 

(4.42f) 



'M M ^/M 3 
We immediately see from eqs.(4.41) that the Fxbb„ are indeed superlunhnal UPWs 
solutions of ME, propagating with speed 1/ cos 77 in the z-direction. That Fxbb„ are 
UPWs is trivial and that they propagate with speed c\ = 1/ cos i] follows because 
Fxbb„ depends only on the combination of variables (z — at) and any derivatives of 
$jbb„ will keep the (z — at) dependence structure. 

Now, the Poynting vector Pxb B n and the energy density uxbb„ for FxBB n are ob- 
tained by considering the real parts of Exbb„ and Bxbb„- We have, 



(PxbbJp = -Re{(ExBBjs}Re{(BxBB n )z} 

(PxbbJb = Re{{ExBB n )p}Re{(B X BB n ) z } 

{PxbbJz = cos?? \\Re{(E X BB n )p}\ 2 + \Re{{E X BB n )e}\ 2 



(4.43a) 
(4.43b) 
(4.43c) 
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u XBBn = (1 + cos 2 77) [\Re{(E X BBjp}\ 2 + \Re{(E XBB J e }\ 2 

+ \Re{(B X BBjz}\ 2 • (4.44) 
The total energy of F x b b„ is then 

£xbb„ = / d0 / dz / P&PU X BB n (4.45) 



Since as 2 — > 00 ExBB n decrease as l/|z — i cos^l 1 / 2 which occurs for the X-branches 
of F X BB n i £xBB n may not be finite. Nevertheless, as in the case of the acoustic X- 
waves, which experiments have shown to travel with v > c s ||, we are quite sure that 
a finite aperture approximation to Fxbl„ (FAAFx bl„) can be launched over a large 
distance. Indeed in |6| computer simulations for the motion of FAAFx bb„ are exhibited 
showing that with an antenna of 20 m of diameter a FAAFx b b„ centered at a frequency 
of 700 GHz propagates with superluminal speed without appreciable distortion up to 
100 Km. See also ^(J. Obviously in this case the total energy of the FAAF X bl„ is 
finite. 

We conclude this Section observing that in general both subluminal and superluminal 
UPW solutions of ME have non-null field invariants and are not transverse waves. In 
particular our solutions have a longitudinal component along the z-axis. This result is 
important because it shows that, contrary to the speculations of Evans p]] ], we do not 
need an electromagnetic theory with a non zero photon-mass, i.e., with F satisfying 
Proca's equation (as proposed also by de Broglie p2fl) in order to have an electro- 
magnetic wave with a longitudinal component. Since Evans presents evidence EJ] of 
the existence of longitudinal magnetic fields in many different physical situations, we 
conclude that the theoretical and experimental study of subluminal and superluminal 
UPW solutions of ME must be continued. 

4.4. The velocity of transport of energy of the UPW solutions of Maxwell 
equations 

Since we found in this paper UPWs solutions of Maxwell equations with speeds < 
v < 00, the following question arises naturally: Which is the velocity of transport of 
the energy of a superluminal UPW (or quasi UPW) solution of ME? 

We can find in many physics textbooks (e.g. |4j|) and in scientific papers |2l|, ^2) the 
following argument. Cons ider an arbitrary solution of ME in vacuum dF = 0. Then if 



F = E + iB (see eq.( 2.77 )) it follows that the Poynting vector and the energy density 
of the field are 

P = ExB, u=^(E 2 + B 2 ). (4.46) 
It is obvious that the following inequality always holds: 

I PI 

v e = — < 1. (4.47) 

u 

Now, the conservation of energy-momentum reads in integral form over a finite vol- 
ume V with boundary S = dV: 
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Eq.(4.48) is interpreted saying that § s dS.P is the field energy flux across the surface 
S = dV, so that P is the flux density — the amount of field energy passing through 
a unit area of the surface in unit time. Now, for plane wave solutions of Maxwell 
equations, 

v £ = 1 (4.49) 

and this result gives origin to the "dogma" that free electromagnetic fields transport 
energy at speed v s = c = 1. 

However v £ < 1 is true even for subluminal and superluminal solutions of ME, as the 



ones discussed in Sections 4.2 and 4.3. The same is true for the superluminal modified 
Bessel beam found by Band EJ in 1987. There he claims that since v e < 1 there is 
no conflict between superluminal solutions of ME and Relativity Theory since what 
Relativity forbids is the propagation of energy with speed greater than c. 

Here we c hallen ge this conclusion. The fact is that as well known P is not uniquely 
defined. Eq.( 4.48| ) continues to hold true if we substitute P i— > P + P' with V.P' = 0. 
But of course we can easily find for subluminal, luminal or superluminal solutions of 
Maxwell equations a P' such that 



\P + P' 



> 1. 



(4.50) 



We arrive at the conclusion that the question of the transport of energy in superluminal 
UPWs solutions of ME is an experimental question. For the acoustic superluminal X- 
wave solution of the HWE (see ^ [|) the energy around the peak area flows together 
with the wave, i.e., with speed c\ — c s j cos 77 while (as we said in Section 3.4) the usual 
theory predicts for the speed of propagation of sound wav es tha t \S\j u < c s , where S 
is the flux of momentum and u is the energy density (eqs. ( |3.8l| ) and (3.82)). Since we 
can see no possibility of the field energy of the superluminal electromagnetic wave to 
travel outside the wave we are confident to state that the velocity of energy transport 
of superluminal electromagnetic waves is superluminal. 



Before ending we give another example to illustrate that eq.(4.47) is devoid of physical 
meaning. Consider a spherical conductor in electrostatic equilibrium with uniform 
superficial charge density (total charge Q) and with a dipole magnetic moment. Then 
we have 

r C 
Q— ; B = —(2 cos 6 r + sin 6 6) 



E 



(4.51) 



and 



Thus 



P = E x B = 



CQ 



2rCQsin9 



2 Q 2 + C 2 (3cos 2 6» + I) 



r 2 

° (1 2 
-^(3 008 



^ 0, for r ^ 0. 



(4.52) 



(4.53) 



Since the fields are static the conservation law eq.(4.3) continues to hold true, as there 
is no motion of charges and for any closed surface containing the spherical conductor 
we have 



dS.P = 0. 



(4.54) 



But nothing is in motion! In view of these results we must investigate whether the 
existence of superluminal UPWs solutions of ME is compatible or not with the Principle 
of Relativity. We analyze this question in detail in the next Section. 
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To end this Section we recall that in Section 2.19 of his book (l4| Stratton presents 
a discussion of the Poynting vector and energy transfer which essentially agrees with 
the view presented above. Indeed he finished that Section with the words: "By this 
standard there is every reason to retain the Poynting-Heaviside viewpoint until a clash 
with new experimental evidence shall call for its revision." 

5. SUPERLUMINAL SOLUTIONS OF MAXWELL EQUATIONS AND THE PRINCI- 
PLE OF RELATIVITY 

In H it was shown that it seems possible with present technology to launch in free space 
finite aperture approximations to the superluminal electromagnetic waves (SEXWs). 
We show in the following that the physical existence of SEXWs implies a breakdown 
of the Principle of Relativity (PR). Since this is a fundamental issue, with implications 
for all branches of theoretical physics, we will examine the problem with great care. In 
Section 5.1 we give a rigorous mathematical definition of the PR and in Section 5.2 we 



present the proof of the above statement. 

5.1. Mathematical formulation of the Principle of Relativity and its physical 
meaning 

In Section || we defined Minkowski spacetime as the triple (M, g, D), where M ~ R 4 , 
g is a Lorentzian metric and D is the Levi-Civita connection of g. Consider now Gm, 
the group of all diffeomorphisms of M, called the manifold mapping group. Let T be a 
geometrical object defined in A C M. The diffeomorphism h G Gm induces a deforming 
mapping h* : T — > /i*T = T such that: 

(i) If / : M D A M, then hj = fo /i" 1 : h[A) -> R. 

(ii) If T G secT^^iA) C secT(M), where T^(A) is the sub-bundle of tensors of 

type (r, s) of the tensor bundle T(M), then 

(h*T)h e (h*ui, . . . , h*u r , h*Xi, . . . , h*X s ) = T e (wi, ...,u r , Xi, ... , X s ) (5.1) 
VXi g T e A, i = 1, . . . , s, Muij G T*A, j = 1, . . . ,r, Ve G A. 

(iii) If D is the Levi-Civita connection and X, Y G sec TM , then 

{h*D h * x KY) he h* f = (D x Y) e f VeeM. (5.2) 



If {f/j. — d/dx^} is a coordinate basis for TA and {# M = dx^} is the corresponding 
dual basis for T*A and if 

T = T£;;£ir- ® . . . ® 6 V > ® ® . . . <g> / Mr , (5.3) 

then 

h*T = [T£:;£ o h-^htff* ® . . . ® /i*^- ® ^/^ ® . . . ® (5.4) 

Suppose now that A and can be covered by the local chart (U, rj) of the maximal 
atlas of M, and A C U, h(A) C [/. Let be the coordinate functions associated with 
(17, rj). The mapping 

2> = 2^ o IT 1 : h(U) R (5.5) 

defines a coordinate transformation (x^) i— > (a; M ) if /i(J7) 3iU Indeed (a; M ) arc 

the coordinate functions associated with the local chart (V, ip) where h(U) C V - and 
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U n V 7^ 4>. Now, since it is well known that under the above conditions h^d/dx^ = 
d/dx" and h*dx>* = dx", eqs.(O), Q and (p|) imply that 



(ft*T) <a / M) (/ie)=T <ie M)(e) J (5.6) 

where T^ M )(e) means the components of T in the chart (x^) at the event e £ M, i.e. 

T (xf .)(e) = ^■■^(^(e)) and where T£"£ r (x> 1 (he)) are the components of T = ft»T 
in the basis {h*d/dx^ = d/dx {/i»d.T M = da; at the point /i(e). Then eg. (5.6) 
reads 

T'Z:t(x»{he)) = T^;(x^e)) (5.7) 



or using eq.(5.£) 

%tt T ^^)) = ( A ~ X )£ ■ ■M'T'£^(x'»(h- 1 a)) (5.8) 

where = dx M /dx a , etc. 

In Section H we already introduced the concept of inertial reference frames I £ sec TU, 
U C M by 

g(I,I) = 1 and £>/ = (5.9) 

A general frame Z satisfies g(Z,Z) = 1, with Z?Z 7^ 0. If a = g(Z, ) £ secT*U, it 
holds 

(.Da)e = a e ® a e + er e + w e + -9 e h e ,e £U C M, (5.10) 

o 

where a = g(A, ), A = DzZ is the acceleration and where w e is the rotation tensor, er e 
is the shear tensor, # e is the expansion and h e — g\ He where 

T e M = [Z e ] 8 [JT e ] • (5.11) 

H e is the rest space of an instantaneous observer at e, i.e. the pair (e, Z e ). Also 
/i e (X,Y) = g e MpX,pY),VX,Y £ T e M and p : T e M -> iJ e . (For the explicit form of 
u>, cr, 0, see |3). From eqs. and ( |57Tot ) we see that an inertial reference frame 
has no acceleration, has no rotation, no shear and no expansion. 

We introduced also in Section || the concept of a (nacs/I). A (nacs/I) is said 
to be in the Lorentz gauge if x^,^ — 0,1,2,3 are the usual Lorentz coordinates and 
I = d/dx° £ secTM. We recall that it is a theorem that putting I = e = d/dx°, 
there exist three other fields £ sec TM, such that g(ej,ej) = —1, i = 1,2,3, and 
e H = d/dx l . 

Now, let (x M ) be Lorentz coordinate functions as above. We say that I £ Gm is a 
Lorentz mapping if and only if 

a>( e ) = A(^(e), (5.12) 

where A£J G is a Lorentz transformation. For abuse of notation we denote the subset 
{£} of Gm such that eq.( 5.12| ) holds true also by C\ c Gjy. 

When (x^) are Lorentz coordinate functions, (x M ) are also Lorentz coordinate func- 
tions. In this case we denote 

e p = d/dx», = a/flx", 7 p = dar", 7 ; = dx'' 1 ; (5.13) 
when t G £1 C G*m we say that £*T is the Lorentz deformed version of T. 
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Let h € Gm- If for a geometrical object T we have 

KT = T, (5.14) 



then h is said to be a symmetry of T and the set of all h G Gm such that eq.( 5.15 ) 



holds is said to be the symmetry group of T. We can immediately verify that for 



S C + C G 



l*g = g, £*D = D, (5.15) 



i.e., the special restricted orthochronous Lorentz group C + is a symmetry group of g 
and D. 

In [fl2| we maintain that a physical theory r is characterized by: 

(i) the theory of a certain "species of structure" in the sense of Boubarki |l5| ; 

(ii) its physical interpretation; 

(iii) its present meaning and present applications. 

We recall that in the mathematical exposition of a given physical theory r, the pos- 
tulates or basic axioms are presented as definitions. Such definitions mean that the 
physical phenomena described by r behave in a certain way. Then, the definitions re- 
quire more motivation than the pure mathematical definitions. We call coordinativc 
definitions the physical definitions, a term introduced by Reichenbach It is nec- 
essary also to make clear that completely convincing and genuine motivations for the 
coordinative definitions cannot be given, since they refer to nature as a whole and to 
the physical theory as a whole. 

The theoretical approach to physics behind (i), (ii) and (iii) above is then to admit 
the mathematical concepts of the "species of structure" defining r as primitives, and 
define coordinatively the observation entities from them. Reichenbach assumes that 
"physical knowledge is characterized by the fact that concepts are not only defined by 
other concepts, but are also coordinated to real objects" . However, in our approach, 
each physical theory, when characterized as a species of structure, contains some implicit 
geometric objects, like some of the reference frame fields defined above, that cannot in 
general be coordinated to real objects. Indeed it would be an absurd to suppose that 
all the infinity of IRF that exist in M must have a material support. 

We define a spacetime theory as a theory of a species of structure such that, if Mod 
r is the class of models of r, then each T £ Mod r contains a substructure called 
spacetime (ST). More precisely, we have 



T = (ST,T 1; ...,T m }, (5.16) 

where ST can be a very general structure fl^| . For what follows we suppose that 
ST = M. = (M,g,D) i.e., that ST is Minkowski spacetime. The Ti,i = 1, . . . ,m are 
(explicit) geometrical objects defined in U C M characterizing the physical fields and 
particle trajectories that cannot be geometrized in T. Here, to be geometrizable means 
to be a metric field or a connection on M or objects derived from these concepts, as e.g., 
the Riemann tensor or the torsion tensor in more general theories. The reference frame 
fields will be called the implicit geometrical objects of r, since they are mathematical 
objects that do not necessarily correspond to properties of a physical system described 
by r. 
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Now, with the Clifford bundle formalism we can formulate in C£(M) all modern phys- 
ical theories (see Section |^) including Einstein's gravitational theory Q . We introduce 
now the Lorentz-Maxwell electrodynamics (LME) in C£(M) as a theory of a species of 
structure. We say that LME has as model 



where (M,g,D) is Minkowski spacetime, {ipi,mi,ei}, i — 1,2,...,N is the set of all 
charged particles, nij and ej being the masses and charges of the particles and ifi : 
M D I — > M being the world lines of the particles characterized by the fact that if 
(fli* £ sec TM is the velocity vector, then Cpi = g(<fi*, ) £ secA 1 (M) C sec C£{M) 
and ipi.ipi — 1. F £ secA 2 (M) C sec C£(M) is the electromagnetic field and J £ 
sec A 1 (M) C sec C£(M) is the current density. The proper axioms of the theory are 



From a mathematical point of view it is a trivial result that tlme has the following 
property: If h £ Gm and if the set of eqs.(5.16) has a solution (F, J, (ipt, rrii, e^)) in 
U C M then {h„F, h*J, (h*ipi, m;, ej)) is also a solution of eqs.(5.16) in h(U). Since the 
result is true for any h £ Gm it is true for £ £ C Gm, i-e. for any Lorentz mapping. 

We must now make it clear that (F, J, {ifi, rrii, ej}) which is a solution of eq.(5.16) 
in U can be obtained only by imposing mathematical boundary conditions which we 
denote by BU. The solution will be realizable in nature if and only if the mathematical 
boundary conditions can be physically realizable. This is indeed a nontrivial point fl^| 
for in particular it says to us that even if (h*F, h^J, {h^ipi, m,-, ej}) can be a solution of 
eqs.(5.16) with mathematical boundary conditions Bh(U), it may happen that Bh{U) 
cannot be physically realizable in nature. The following statement, denoted PRi, is 
usually presented (l2) as the Principle of (Special) Relativity in active form: 

PR\\ Let £ E C\_ C Gm- If for a physical theory r we have T S Modr, where T = 
(M, g, D, Ti, . . . , T m ) is a possible physical phenomenon, then £*T = (M, g, D, l*Ti, 
. . . , Z*T m ) is also a possible physical phenomenon. 

It is clear that hidden in PR\ is the assumption that the boundary conditions that 
determine £*T are physically realizable. Before we continue we introduce the statement 
denoted PR2 known as the Principle of (Special) Relativity in passive form Q. 

PR2: All incrtial reference frames arc physically equivalent or indistinguishable. 

We now give a precise mathematical meaning to the above statement. 

Let r be a spacetime and let ST = (M, g, D) be a substructure of Mod r representing 
spacetime. Let I £ secTU and I' £ sec TV, U, V C M be two inertial reference frames. 
Let (U, 77) and (V, tp) be two Lorentz charts of the maximal atlas of M that are naturally 
adapted respectively to / and For (x 11 ) and (x M ) the coordinate functions associated 
with ([/, rj) and (V, ip) we have / = d/dx°, I' = d/dx °. 

Definition. Two inertial reference frames / and i 7 as above are said to be physically 
equivalent according to r if and only if the following conditions are satisfied: 



(i) Gm D £+ 3 £ :U -> £{U) CV, x'^ = x^o £~ x =>- V = £J 
When T £ Mod r, T = (M, g, D, Ti, . . . T m ), is such that g and D are defined over 
all M and T. ( £ secC£(U) C secC^(Af), i = l,...,m, calling o = (g, D, Ti, . . . T m ), 



Tlmb = (M, g, D, F, J, {<fii,mi, ej 



}) 



(5.17) 



dF = J , 
m.iD^pi = e t Cpi ■ F . 



(5.18) 
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o solves a set of differential equations in f](U) C M 4 with a given set of boundary 
conditions denoted b°^ x "\ which we write as 

£>f x , ) (o<x.>)e=0 ; 6°<^> - eel! (5.19) 

and we must have: 

(ii) If T e Mod T»f,Te Mod r, then necessarily 

4T = (M l5 ,A4Ti,..i«T ro ) (5.20) 

is defined in £([/) C V and calling = {g, D, £*Ti, . . . , £*T m } we must have 

The system of differential equations (5.19) must have the same functional form as the 

system of differential equations (5.17) and b e *°( x ) must be relative to (x M ) the same 

as b ^^ is relative to (x^) and if b ^^ is physically realizable then b e * o! < x ) must also 
be physically realizable. We say under these conditions that I <~ /' and that £*o is the 
Lorentz version of the phenomenon described by o. 

Since in the above definition £*Y = (M, g, D,l*T\, . . . ,£*T m ), it follows that when 
/ ~ I', then £ 3t g — g, £*D — D (as we already know) and this means that the spacetime 
structure does not give a preferred status to I or I' according to r. 

5.2. Proof that the existence of SEXWs implies a breakdown of PRi and 
PR 2 

We are now able to prove the statement presented in the beginning of this Section, that 
the existence of SEXWs implies a breakdown of the Principle of Relativity in both its 
active (PRi) and passive (PR2) versions. 

Let £ e C\ C G M and let F, F E secA 2 (M) C aecC£(M), F = £*F. Let F = 
LF = RFR- 1 , where F e = (l/2)F fny (x^(£- 1 e))^-f l/ and where R e sec Spin+(1,3) C 
sec C£(M) is a Lorentz mapping, such that 7'^ = Rj^R^ 1 = A^7 Q , A£ e C\ and let 
(x^) and (x M ) be Lorentz coordinate functions as before such that 7^ = dx^, 7 ^ = di 11 
and i^i^o £~ x . We write 

F e = l -F,Ax»{e))n v (5.22a) 
Fe = |^(i'"(e)) 7 V (5.22b) 
F e = ^(^(eJJYV (5.23a) 

F e = ^(;"(e)) 7 V. (5.23b) 
From (5.22a) and (5.22b) we get that 



F^(x»(e)) = (A- 1 r a (A-%F^(x»(e)). 



(5-24) 
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From (5.23a) and (5.23b) we also get 



(5.25) 



Now, suppose that F is a superluminal solution of Maxwell equation, in particular 
a SEXW as discussed in Section 3. Suppose that F has been produced in the inertial 
frame I with (x^) as (nacs/7), with the physical device described in Section 3. F is 
then traveling with speed c\ = 1/ cos r\ in the negative z-direction and being generated 
in the plane z = 0, will travel to the future in spacetime, according to the observers in 
I. Now, there exists I € C\_ such that £*F = F = RFR^ 1 will be a solution of Maxwell 
equations and such that if the velocity 1-form of F is vp = {c\ — 1) -1 / 2 (1, 0, 0, — ci), 
then the velocity 1-form of F is Vp = (cj 2 — l) _1 / 2 (— 1, 0, 0, — c[), with c[ > 1, i.e, v-p is 
pointing to the past. As its is well known F carries negative energy density according 
to the observers in the I frame. 

We then arrive at the conclusion that to assume the validity of PRi is to assume the 
physical possibility of sending to the past waves carrying negative energy. This seems 
to the authors an impossible task, and the reason is that there do not exist physically 
realizable boundary conditions that permit the observers in I to launch F in spacetime 
and such that it travels to its own past. 

We now show that there is also a breakdown of PR2, i.e., that it is not true that 
all inertial frames are physically equivalent. Suppose we have two inertial frames I and 
V as above i.e, I = d/dx Q , I' = d/dx °. Suppose that F is a SEXW which can be 
launched in I with velocity 1-form as above and suppose F is a SEXW built in V at 
the plane z' = and with velocity 1-form relative to (x M ) given by v-p — v and 



If F and F are related as above we see that F, which has positive energy and is traveling 
to the future according to can be sent to the past of the observers at rest in the / 
frame. Obviously this is impossible and we conclude that F is not a physically realizable 
phenomenon in nature. It cannot be realized in V but F can be realized in /. It follows 
that Pi?2 does not hold. 

If the set of inertial reference frames are not equivalent then there must exist a 
fundamental reference frame. Let I £ sec TM be the fundamental frame. If F is 
moving with speed V relative to I, i.e. 



then, if observers in V are equipped with a generator of SEXWs and if they prepare 
their apparatus in order to send SEXWs with different velocity 1-forms in all possible 
directions in spacetime, they will find a particular velocity 1-form in a given spacetime 
direction in which the device stops working. A simple calculation yields then, for the 
observers in /', the value of VI 

In Recami argued that the Principle of Relativity continues to hold true even 
if superluminal phenomena exist in nature. In this theory of tachyons there exists, of 
course, a situation completely analogous to the one described above (called the Tolman- 
Regge paradox), and according to Recami's view PR2 is valid because V must interpret 
F a being an anti-SEXW carrying positive energy and going into the future according to 




(5.26) 



/' = 



15 V d 



(5.27) 



VI - V 2 dt VI - V 2 dz ' 
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him. In his theory of tachyons Recami was able to show that the dynamics of tachyons 
implies that no detector at rest in / can detect a tachyon (the same would be valid for 



a SEXW like F) sent by I' with velocity 1-form given by eq.(5.26). Thus he claimed 
that Pi?2 is true. At first sight the argument seems good, but it is at least incomplete. 
Indeed, a detector in / does not need to be at rest in /. We can imagine a detector in 
periodic motion in / which can absorb the F wave generated by /' if this was indeed 
possible. It is enough for the detector to have relative to I the speed V of the /' frame 
in the appropriate direction at the moment of absorption. This simple argument shows 
that there is no salvation for PR2 (and for PRi) if superluminal phenomena exist in 
nature. Our argumentation is endorsed by Barashenkov and Yur'iev |48[ . 

The attentive reader at this point probably has the following question in his/her 
mind: How could the authors start with Minkowski spacetime, with equations carrying 
the Lorentz symmetry and yet arrive at the conclusion that PR\ and PR2 do not 
hold? The reason is that the Lorentzian structure of (M,g,D) can be seen to exist 
directly from the Newtonian spacetime structure as proved in [^9). In this paper, 
Rodrigues and collaborators show that even if £ \ is not a symmetry group of Newtonian 
dynamics it is a symmetry group of the only possible coherent formulation of Lorentz- 
Maxwell electrodynamic theory compatible with experimental results that is possible 
to formulate in the Newtonian spacetime J^| 

We finish calling to the reader's attention that there are some experiments reported in 
the literature which suggest also a breakdown of PR2 for the roto-translational motion 
of solid bodies. A discussion and references can be found in Jl^|. A coherent spacetime 
model which can accommodate superluminal phenomena has been recently proposed 
by Matolcsi and Rodrigues fsojj . 

6. CONCLUSIONS 

In this paper we presented a unified theory showing that the homogeneous wave equa- 
tion, the Klein-Gordon equation, Maxwell equations and the Dirac and Weyl equations 
have solutions with the form of undistorted progressive waves (UPWs) of arbitrary 
speeds < v < 00. We exhibit also some subluminal and superluminal solutions of 
Maxwell equations. We showed that subluminal solution can in principle be used to 
model purely electromagnetic particles. 

The possible existence of superluminal electromagnetic waves implies in a breakdown 
of the Principle of Relativity. It is important to recall here that exact Lorentz symmetry 
can be preserved in an abstract mathematical level through the ingenious construction 
of Santilli's isominkowskian spaces (see [51-55]). Santilli's theory is important, e.g. for 
situations involving the hadronic medium, where superluminal velocities can occur. We 
observe that besides the fundamental theoretical implications, the practical implications 
of the existence of UPWs solutions of the main field equations of theoretical physics 
(and their finite aperture realizations) are very important. This practical importance 
ranges from applications in ultrasound medical imaging to the project of electromag- 
netic bullets and new communication devices |5(|. Also, we would like to conjecture 
that the existence of subluminal and superluminal solutions of the Weyl equation may 
be important to solve some of the mysteries associated with neutrinos. Indeed, if neu- 



trinos can be produced in subluminal or superluminal modes (see ]57j, 58 for some 
experimental evidence concerning superluminal neutrinos) they can eventually escape 
detection on earth after leaving the sun. Moreover, for neutrinos in a subluminal or 



2 We recall that Maxwell equations have, as is well known, many symmetry groups besides C,. 
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superluminal mode it would be possible to define a kind of "effective mass" . Recently 
some cosmological evidences that neutrinos have a non-vanishing mass have been dis- 
cussed, e.g., by Primack et al p9| . One such "effective mass" could be responsible for 
those cosmological evidences, and in such a way we can still have a left-handed neutrino 
since it would satisfy the Weyl equation. We shall discuss more this issue in another 
publication. 
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